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Riassunto. — Sono state fatte misure in alta frequenza (w = 4,39 : 10°) 
della costante dielettrica dell’acqua quando in essa siano stati disciolti 
degli elettroliti. Si & trovato che la costante dielettrica dell’acqua dimi- 
nuisce e poi aumenta al crescere della concentrazione c, riscontrandosi 
il minimo in corrispondenza di e dell’ordine di 10°. 


Le misure della costante dielettrica dell’acqua in soluzioni elettrolitiche 


non hanno ancora raggiunto una conclusione sicura; esse sono state attuate 


da vari sperimentatori, coi metodi più diversi, e con risultati non molto in 


accordo fra loro. In questa Nota noi 
vogliamo riferire di aleune misure at- 
tuate, su tale argomento, con un me- 
todo che si presenta semplice, e che 
quindi permette la raccolta, in modo 
rapido, di parecchi dati relativi a 
varie soluzioni. Esso è basato sul fe- 
nomeno di risonanza in un circuito 
in cui c’é una resistenza R ed una 
capacità C da misurare; la fig. 1 lo 
riproduce schematicamente. Z è la 
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bobina che permette un collegamento induttivo con un oscillatore di fre- 
queuza » — 7 MHz non rappresentato in figura; la cella contenente la soluzione 
elettrolitica € rappresentata schemati- 
camente da una capacita €; e da una 
resistenza R posta in parallelo; in 
queste misure puo parlarsi esplicita- 
mente di capacità della ceila, perchè 
sono stati adoperati due elettrodi 
estesi e molto vieini (la capacità cal- 
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colata approssimativamente dai dati 
sof geometrici & di circa 1,6 pF). L’autoin- 
duzione L, à accoppiata ad un circuito 

in cui ¢ un diodo rettificatore che per- 
mette di misurare una intensità di cor- 

rente proporzionale a quella che circola 

nel cireuito di fig. 1. Facendo ora 
variare C, si raggiungerà il massimo 
della corrente del cireuito nelle condi- 

D zioni di risonanza. Evidentemente per 
AG 25 aA ogni soluzione posta nella cella si ha 
Fig. 2. uno speciale valore della capacita C, e 

della resistenza R e quindi una deter- 

minata posizione della capacità C ed un determinato valore della corrente, 
nel massimo di risonanza. Ora, per fare la misura, si pone nella cella €, una 
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serie di soluzioni di acetone in acqua con costante dielettrica A nota e, per 
ogni soluzione, una serie di resistenze, derivate agli elettrodi, di valore com- 


preso fra 50 e 1000 ohm, in modo da 
conoscere per ogni Singolo valore della 
capacità ©, e della resistenza R la po- 
sizione D, del condensatore C nel mas- 
simo della corrente (condizione di ri- 
sonanza) ed il valore D,, Jetto al gal- 
vanometro, proporzionale al valore 
della intensitä di corrente, in tali con- 
dizioni di risonanza; si ottengono, ri- 
portando i valori di D, e di D,, due 
famiglie di curve: le une sono quelle 
che si ottengono mantenendo costante 
la resistenza R e facendo variare la 
eapacitä C,; le altre mantenendo in- 
vece costante la capacità C, e facendo 
variare la resistenza R (nelle figure 
sono indicate con curva tratteggiata 
le isoohmiche, e con curva continua 
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le isofarad). Costruito in tal modo questo reticolato di curve, si pud porre 
nella cella una serie di soluzioni saline di varie concentrazioni, in modo da 
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avere soluzioni a cui competono 
certi valori di R, e se per la pre- 
senza del sale la costanza dielet- 
trica dell’acqua non varia, i 
punti figurativi delle soluzioni 
con diversa concentrazione, e 
cioé con diversa R, si devono 
trovare su di una linea con an- 
damento compreso fra due suc- 
cessive curve isofarad; in caso di- 
verso, e cioè se c’é variazione 
di K, la curva deve avere un an- 
damento non parallelo alle iso- 
farad di taratura. 

Le misure compiute in un pri- 
mo ciclo sperimentale sono rela- 


tive alle soluzioni di KCl (1-1), Na,SO, (1-2), KH,PO, (1-3), K,FeCN, (1-4) e 
di un acido debole, CH,COOH, le cui concentrazioni, come ordine di grandezza, 
vanno da 1/100 di mole/litro ad 1/1000 di mole/litro. Le fig. 2, 3, 4, 5, 6 ripor- 
tano i risultati ottenuti; da esse risulta che le soluzioni elettrolitiche presen- 
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tano effettivamente una variazione di costante dielettrica al variare della con- 
centrazione dell’elettrolito. Dalla fig. 2, per esempio, nella quale, una volta 
per tutte, sono segnati i valori di X e di R nella famiglia delle curve di tara- 
tura, e che si riferisce alle soluzioni di KCl, & da notare, più particolarmente, 
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che le soluzioni più diluite, e cioè quelle a cui compete una R più grande, 
danno luogo a punti figurativi che appartengono a curve a cui corrisponde 
un certo valore di K, mentre i punti rappresentativi di soluzioni più concen- 
trate (e cioè per R più piccolo), vanno a cadere su di una curva a cui corri- 
sponde un valore di K inferiore al precedente, per riappartenere ad una iso- 
farad di valore più vicino al primo per soluzioni di ancora più alta concen- 
trazione. Dai punti figurativi si pud, con una certa approssimazione, ricavare 
l’andamento del valore di K al variare della concentrazione. Tale andamento 
è riportato in fig. 7 per il K,FeCN,. Misure quindi eondotte con metodo com- 
pletamente diverso da quello usato da R. Ftrru (!) hanno portato ad eguale 
risultato: una diminuzione e poi un aumento della costante A al crescere di e. 


(4) R. FÜRTH: Phis. Zeits., 25, 676 (1924). 


SUMMARY 


Measures have been made at high frequency (@ = 4.39 107) of the water dielectric 
constant when some electrolyte has been dissolved in it. It has been found that 
increasing the concentration c, the dielectric constant of the water diminishes and 
afterwards increases and the minimum has been ascertained to be when e is in the 
range of 10°, 
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Summary. — The experimental results on the composition of cosmic 
radiation at the top of the atmosphere are used to develop a pheno- 
menological theory of all processes which occur in the cosmic radiation 
while it traverses the atmosphere. We have computed the distribution 
of the fast nucleonic component (protonic and neutronic), the distri- 
bution of the mesonie component and its energy spectrum and positive 
excess at various altitudes, and the distribution of the electronic com- 
ponent. The latitude effects due to the earth’s magnetic field are studied 
accurately. The theoretical. results obtained are discussed, in the light 
of the most recent experimental data. 


Introduction. 


The primary component of cosmic radiation, which in traversing the 
atmosphere gives rise to most of the processes subject to experimental inve- 
stigation, was thought (1) (?), until recently, to be formed mainly of protons. 


(*) This paper is an organic exposition of the results obtained by one of the authors 
(P. CALDIROLA) and his coworkers during the last two years. Some of these results 
have already been published in preliminary papers: 
P. CALDIROLA: Atti del Congresso Internazionale di Como, in Suppl. Nuovo Oimento, 
6, 511 (1949). 
P. CALDIROLA: Nuovo Cimento, 6, 565 (1949) (referred to as I). 
P. CALDIROLA and A. LOINGER: Nuovo Cimento, 7, 161 (1950) (II). 
P. CALDIROLA and G. Zin: Nuovo Cimento, 7, 575 (1950) (IIT). 
P. CALDIROLA and P. GULMANELLI: Nuovo Cimento, 8, 229 (1951) (IV). 
P. CALDIROLA, R. Frescur and P. GULMANELLI: Nuovo Cimento, 8, 508 (1951) (V). 


: () T. H. Jonnson: Rev. Mod. Phys., 11, 208 (1939). 
(2) M. Schein, W. P. Jesse and E. O. Worran: Phys. Rev., 59, 615, 930 (1941). 
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More recent experiments (?) have proved that protons represent about 70— 75% 
of the primary radiation, while the remaining part is formed essentially of 
a-particles and of a smaller number of heavier nuclei. 

The main problem in cosmic ray theory is to predict the existence and the 
behaviour in the atmosphere of the various components which arise from the 
primary protons and nuclei; and several attempts have already been made 
in this direction. Here we shall tackle the problem with the simplest hypo- 
thesis by studying in detail the following questions: a) the variation of the 
nucleonie component through the atmosphere; b) the origin and behaviour 
of the mesonic component, with emphasis on the existence of a positive excess 
varying with energy, altitude and latitude; and ec) the origin and behaviour 
of the electronic component. We shall also study the latitude effect on the 
nucleonic and mesonic components, since this effect is connected with the 
process of generation of mesons. 

Two different theories about meson production have been worked out 
some time ago: the multiple production theory due to WATAGHIN (4), HEISEN- 
BERG (°) and OPPENHEIMER LEWIS and WOUTHUISEN (°); and the plural pro- 
duction theory due to HEITLER and JANOSSY (7). Recent experiments have 
lead (5) physicists to believe that meson production has an intermediate cha- 
racter: in the collision between two nucleons the production is practically 
single up to a certain energy and then becomes multiple with increasing mul- 
tiplicity for higher energies. A more satisfactory theory of meson production 
may then be the one recently developed by FERMI (°). 

The scheme that we have adopted since 1949 is essentially phenomeno- 
logical, though qualitatively based on the meson theory of nuclear forces. 
Indeed this theory in its «symmetrical » form justifies qualitatively the main 


(3) P. FREIER, E. L. LONGFREN, E. P. Ney and F. OPPENHEIMER: Phys. Rev., 
74, 1818 (1948); H. L. Brapr and B. PETERS: Phys: Rev., 74, 1838 (1948). 

(4) G. WATAGHIN: Phys. Rev., 74, 975 (1948). 

(5) W. HEISENBERG: Nature, 164, 65 (1949); Alti del Congr. Int. di Como in Suppl. 
Nuovo Cimento, 6, 491 (1949). 

(5) H. W. Lewis, J. R. OPPENHEIMER and S. A. WOUTHUISEN: Phys. Rev., 73, 
127 (1948). 

(7) W. HEITLER and L. Janossy: Proc. Phys. Soc., A 62, 374 (1949); A 62, 669 
(1949). The theoretical foundations of the hypothesis assumed here are exposed in 
the following papers: J. HAMILTON, W. Heitrpr and H. W. PENG: Phys. Rev., 64, 
78 (1943); W. Hernrer and P. WALSH: Rev, Mod. Phys., 17, 252 (1945); L. JANOSSY: 
Phys. Rev., 64, 345 (1945). Recently the scheme of HEITLER-JAanossy has been further 
developed by MESSEL, particularly concerning low energies: Communications of the 
Dublin Institute for Advanced Studies, Series A. No. 7 (1951). 

(*) See papers presented at the Como Conference (Sept. 1949) and relative discussions: 
Atti del Congr. Int. di Como in Suppl. Nuovo Cimento, 6, 310 (1949). 

(?) EK. FERMI:: Progress in theor. Phys., 5, 570 (1950); Phys. Rev., 81, 683 (1951). 
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features of the processes arising from the interaction of nucleons and its fre- 
quent failure to yield quantitative results is probably connected (2°) with the 
inapplicability of current quantum methods of calculation to interactions 
between particles nearer than the elementary length A, = hju,c (u, —.mass 
of the meson coupled to the nucleons identified with the x-meson). We assume 
that nuclear interactions are qualitatively of the symmetrical meson type and 
we derive by semiempirical considerations hypotheses apt to specify quan- 
titatively the laws of interaction. In the collision of two nucleons of suffi- 
ciently high energy one or more mesons are created: at the beginning no 
hypothesis is made on the number of mesons, on their distribution with energy 
and among the three types z+, x, z°. Indeed, to derive the distribution law 
of the nucleonie component in the atmosphere, it is sufficient to specify how 
the kinetic energy of the incoming nucleon is distributed among the nucleon 
knocked-on, the scattered one and the mesons created. The comparison of 
the results with experimental data leads one then to specify the mechanism 
of creation of mesons and so to adopt hypotheses which allow to describe the 
process analitically and to calculate the distribution law of mesons in the 
atmosphere. The distribution of the electronic component is calculated with 
the theory of electromagnetic showers. 


1. — Analysis of the nueleonie component. 


11. — Interaction between nucleons. 


We shall make rather simple assumptions on the energy distribution in 
the collision of two nucleons with the production of z-mesons in order to 
obtain expressions amenable to numerical calculation. Our hypotheses are 
taken from a paper by HOcKER (1). The incident nucleon with kinetic energy E 
gives up in a single collision a fraction e*H, (e* constant) of its energy to the 
mesons created, while the remaining energy (1 — e*)H, = ef, is given to the 
two nucleons; owing to the high energy of the nucleons involved, we shall 
assume that the scattered nucleon can assume with equal probability any 
amount H of energy between 0 and eH, so that the nucleon knocked-on will 
assume the energy eH, — LH. 

Energetically the creation of mesons in a nucleon-nucleon collision can 
occur if the energy of the incident nucleon is at least equal to the threshold 
energy uw, Obviously, if the energy is very near to the threshold the cross 
0500) W. HEISENBERG: Zeits. f. Phys., 101, 533 (1936). 

(1) K. H. Höcker: Zeits. f. Phys., 124, 351 (1948). 
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section is very small, as the artificial production of mesons proves (2): the 
cross section will increase with increasing energy toward the geometrie section 
as one observes in cosmie ray phenomena. For sake of simplicity we shall 
assume that the cross section o is « effectively » different from zero only for 
values of the kinetie energy of the incident nucleon higher than or equal to 
a critical value Z,; and that for these values of the energy o is effectively 
constant. We shall assume also that for energies > E, the nucleons lose 
energy only in collision processes with other nucleons: indeed, since E, — 
— +3 GeV (GeV = 10° eV: gigaelectronvolt), the ionization loss for protons 
is small and the collision of nucleons with nuclei leading to the usual disinte- 
grations has a much smaller cross section. 


1°2, — Distribution of fast nucleons in the atmosphere. 


We need to consider only the distribution of nucleons with kinetic energy # 
greater than or equal to the critical value Z,: we shall call these nucleons 
«fast » (8). If the cut-off energy of the earth’s field E,< E,, there is no 
latitude effect on the production of mesons and on the subsequent processes. 
Here we shall study the global distribution of nucleons subject to the hy- 
pothesis HE, < #,: such a distribution can be calculated without further spe- 
eifications on meson production. 

H(E, 1) is (*) the total number of nucleons having an energy >E at a 
depth / in the atmosphere, measured in units characteristic of the nucleon- 
nucleon collision (5); h(#, 1) is the corresponding spectral density. So 

eo 
(1) H(E, 1) = | h(E',D dE:  h(E, 1) = 


E 


_ 2H, 1) 
dE 


(2) E. GARDNER and C. M. G. Larres: Science, 107, 270 (1948); J. BURFENING, 
E. GARDNER and C. M. G. LATTES: Phys. Rev., 75, 382 (1949); V. Z. PETERSON: Phys. 
Rev., 79, 407; 80, 136 (1950); W. F. CARTWRIGHT, C. RICHMAN, M. A. Wırcox, and 
M. H. WHITEHEAD: Phys. Rev., 78, 823; 79, 198 (1950); 81, 652 (1951). 

(3) B. FERRETTI: Nuovo Cimento, 6, 379 (1949) has divided the nucleonie com- 
ponent into three classes: class A, including nucleons with an energy higher than 
600 MeV; class B, including nucleons with an energy between 200 and 600 MeV; and 
class C including nucleons with an energy less than about 200 MeV. The fast nucleonic 
component considered by us is formed by the most energetic fraction of the nucleons 
belonging to class A. This component is the most efficient source of mesons. 

(4) The calculations of this section follow those of HÖCKER. 

(5) If we indicate with h the vertical height in centimeters, with o the nucleon- 
nucleon cross section in cm?, with A = A, exp (—yh) (y = 1/804000) the number of 
nucleons contained in one cm? of atmosphere, we have di = — A(h)dh: where from 
1 = (A55/y) exp (— yh). 
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The conservation equation reads: 
(2) H(E, 1+ dt) = H(E,1) + dH, (#, 1) — dH.(E, 1) , 


if dH, and dH. are the numbers of nucleons which in the depth di have 
increased their energy above E, or, decreased it below E, respectively. With 
the assumptions of 1°1 and using the one-dimensional approximation (the 
scattered particle and the one knocked-on have the same direction as the 
incident particle) we have: 


Ele = 
dH.(E, 1) = | AB’: h(E', 1) dt + | dE'-h(E", Da, 
E Ele 
> 
aH-(E, l) = | dE’: h(B’, 1) A Het 
va j eH 


Then equation (2) yields, by partial integration, the following integro-diffe- 
rential equation: 


Eu aa nee 
(3) Zn = H(B,D—° | dE HE, |) zn: 
Ele 


The primary component at the top of the atmosphere is formed of protons 
with an energy spectrum P(E, 0) = P,: E- and of heavy nuclei with a total 
number of nucleons: 


H*(E, 0) = P*(H, 0) + N*(E, 0), 


where P* and N* are approximately equal. For the sake of simplicity we 
shall assume for these nucleons the same energy spectrum as for the primary 
protons: 

HERO) = HB. 


Since heavier nuclei are less effective than protons to produce mesons by 
interaetion with the nuclei of the air (also because they lose more energy by 
ionization), we replace the primary nuclear component with an effective com- 
ponent of free nucleons, by putting: 


(4) H(E,0) = P(E, 0) + AH*(E, 0) = (P, + AH*)E- = HE, 
where A<1. The solution of (3), with the condition (4), is: 
DES 
(5) H(E, 1) = H,E" exp (— al) e nl i) ; 
. \ It, 


This is the distribution law of fast nucleons in the atmosphere. The agree- 
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ment with experiment is good if s is taken between 1.5 and 2 (we shall choose 
1.8 as most authors do) and « is taken to give a mean free path for absorption L 
between 110-130 g/cm? (1%). The distribution law of nucleons in the atmo- 
sphere can as well be written as: 


(5)  H(E,x) = H,E- exp (— w/a») exp (w/a,) = H,E- exp (—a/L), 


and the value of x, can be determined from the interaction of nucleons with 
air nuclei (7). Considering also the experimental results (18) in materials of 
atomic weight comparable to that of the air nuclei, one finds that a) is between 
60 and 80 g/cem?; the cross section is then slightly smaller than the geometric 
one («transparency » of nuclear matter) (7%). 

The exact value of x, cannot be found by direct measurements, since in 
our theory air is schematized as an homogeneous distribution of nucleons. 
The values mentioned for x, and Z give, through the identity of equations (5) 
and (5') 


ze iy 
Ser 


a value between 0.5 and 0.9 for e: a more exact estimate of e is provided by 
latitude effects (see n. 6). In our calculations we have assumed Z = 125 g/cm? 
and e — 3/4. In the collision between a nucleon of kinetic energy Æ and a 
nucleon of the atmosphere, a quarter of Æ is given to the mesons produced. 


2. — Mechanism of meson production. 


Let us consider the multiplicity k of meson production in a collision between 
two nucleons: k will be a function of the energy E of the incident nucleon. 
The mesonic component of cosmic rays shows a rather definite latitude effect 
at 9 = 45°: the corresponding cut-off of the proton spectrum equals 3 GeV (°°). 


(6) G. BERRNADINI, G. CORTINI and A. MANFREDINI: Phys. Rev., 74, 1878 (1948); 
G. WATAGHIN: Phys. Rev., 71, 453 (1947): I. TınLoT: Phys. Rev., 73, 1456 (1948); 
E. P. GEORGE and A. P. JAsSon: Nature, 161, 248 (1948); R. MAZE, A. FRÉON, J. DAUDIN 
and P. AUGER: Rev. Mod. Phys., 21, 14 (1949). 

(7) G. WATAGHIN: Phys. Rev., 75, 693 (1949). 

(8) G. CocconI: Phys. Rev., 76, 984 (1949). 

(%) The value assumed for the cross section of nucleon-nucleon interaction deter- 
mines the value of the thickness of atmospheric matter corresponding to ? = 1. From 
the assumed value of L, it follows that ! = 1 is equivalent to about 72 g/em?. Con- 
sequently the whole atmospheric depth equals 14.4. 

(22) See L. JANossy: Cosmic Rays (Oxford, 1948), page 306. 
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This suggests that nucleons with an energy of few GeV are already capable 
of producing mesons: actually most of the mesons must be generated from 
primary nucleons with energy between 10° and a few times 101 eV since the 
primary spectrum decreases rapidly with energy. In this limited field of 
energy we may assume for k a constant value: a nucleon of energy E creates 
in a single collision & mesons 7 with mean energy (intrinsic + kinetic) 
Dr EN (En) 

OCCHIALINI, POWELL and coworkers (??) have shown that the 7-mesons 
transform into p-mesons with a total mean energy of about 4/5 of that 
of the generating «-mesons: if one takes Æ between 1 and 10 GeV, the mean 
energy of the mesons generated by the primary nucleons will be between 
0.2/k GeV and 2/k GeV. Therefore in this interval of energy k cannot be 
much greater than 1: otherwise mesons derived from protons of few GeV 
would be almost absent from the meson spectrum at low altitude in contrast 
with experimental evidence. Several theoretical and experimental arguments 
had already been presented to suggest that, for £ 2 10° = 10 eV, k is about 
unity (2%): meson production would be predominantly due to plural processes. 

We shall assume k — 1: this must be considered, however, as an hypothesis 
to be checked by its consequences on the behaviour of the nucleonic and 
mesonic components in the atmosphere. Furthermore the assumption will 
not be valid at high energies. According to the symmetrical theory of meson 
fields we shall assume that in the nucleon-nucleon collision the following cases 
ean occur. 

a) proton-proton collision; two possibilities exist: 


| PEN + rt (probability y) 


1er ee 
| 2 + P + ( »  1—» 


b) neutron-neutron collision; two possibilities exist: 
PENIT (probability y) 
Ay Se Iai 


N + 7° ( » 1 =) 


(CS) MIE experimental results (U. CAMERINI, P. H. FOWLER, Ww. ©. Lock and H. 
Murrunap: Phyl. Mag., 41, 413 (1950)) on a power speetrum in the production of mesons 
do not eontradiet this assumption. Indeed such a speetrum implies that the mean 
energy of a meson is proportional to the energy of the incident nucleon. 

(22) C. M. G. LATTES, Gr. P. $S. Occnrarını and C. F. POWELL: Nature, 160, 413 
and 486 (1947). The fraction of energy transferred to the meson is the same, if the 
meson is in motion or at rest: this can be proved by a Lorentz transformation. 

(23) Various arguments supporting a prevalent produetion of mesons through plural 
processes have been presented at the Como Conference (Sept. 1949) by L. JANOSSY 


and B. FERRETTI. 
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e) neutron-proton collision; three possibilities exist: 


& [ Pepe (probability v) 
eae ales v) 
ae | pee no ( » 1— Br). 


3. — Distribution of fast protons and neutrons in the atmosphere. 


311. — Solution of the system of integro-differential equations. 


The distribution functions P(E, 1), p(E, 1!) for the fast protons and N(#, l), 
n(£,1) for the fast neutrons are calculated using the scheme of reactions. 
given above and the procedure of HOCKER (12), properly extended. 

Let 7 be the percentage of protons among all nucleons composing the nuclei 
in the atmosphere, then 


Ble CO 
aP. Vf ' AI nd 
un | PUB Dae + (Lay | pe’, D de 
E Ble 
+ m+ ma (pe, 0 À am 
Kl a a 1] pe, b Arie , 
Ele 
a Wr eE' E al 
ar en Gap] | ote, Tr dE 
Ble 
[ eE'— E 
+ [A —n)y + 9] } ne, D dE 
Ele ; 
and 
Ele co 
an — ; 
edie", [ UL", DAE + m | n(£", 1) dE" + 
E Ele 
A E 
ir ee a EI OLB 
: eh 
Bis 
an. : Ar eE'— E 1 
ne = [d= y) + op) | OB) = a+ 
Ble 


CESR 
+ by + | pee D 2 a’. 


Ele 
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The conservation equation yields by partial integration the following system 


of equations: 


| arm 4 Tps da 
Ele 
ai 
a NED pe à 
S Ele 
(6) ie 
IN (E, I) | AC 9.5 
En ne NE | ; NE, De — 
SS = ME, DU + Ami | XD ge 
Ele 
E | dE 
— [ny — (A — PE) 
| mol: | PIE) on 
Ele 


Obviously by summing these two equations one gets equation (3). 
To apply equations (6) to the distribution of the protonie and neutronic 


components in the atmosphere, we may assume 7 = 1/2; then: 
= eo 
dP(E, |) af E 4 fee Lig ary 
eee ee P(E, I= | PIE), IRAN (DR a Oe 
dl P(E l) ae ( ) iB is ( In 2 
Ele Ele 
1 
where a=1-7 5 (1 — y) 
6") = Rothe 
dN(E, |) € APE ES 0 | Fe vie 
EU WN (EN) —— | NEN di Al PE, bin, dE 
= N(E, D AE Dee 5 | PUD ga 
Ble Ble 


1 
where b = 5 (yap L). 
The solutions of (6’), which satisfy the initial conditions 


PE, 0) = (Py AnS/2)H +,  N(E, 0) = (AHS /2)E , 


are: 
| E- % 4 = Het) 
Pt) ="; [Hoexp (—al) + Pop EAN | ae 
a Ai wt wi emi 
N(E, l) = 5 [Ho exp Ge al) — PS exp — pt] | p wll tree: ; 4 


14 pP! 


So the spectral distribution functions are: 
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il 
= 5 SEA, exp (— al) + P, exp (— Bl)] 


1 
= 5 SE TA, exp (— al) — P, exp (— fl)] 


b* P(E, |) 
| p(k, l) == > JE 
(8) | n(E, l) = IE 
H(E, | 
| ME, 1) = \ a D) = H,sE-* exp (— al) . 


There are no direct experiments which allow to determine the probabili- 
L p | 


~ 


2c 


Nfem?.s.sterad 


ties y and v. We have as- 
sumed y = 1/2 and y = 1/3 
thus considering the diffe- 
rent possibilities for each 
type of interaction equally 
probable for sake of simpli- 
city and owing to the sym- 
metry of nuclear forces (24): 
also some later work of 
QUERCIA, RISPOLI and SCIU- 
TI (?*) on the dependence of 
the positive excess of the 
mesonic component on alti- 
tude and direction upholds. 
the values chosen. These 
values yield « = 0.57 and. 


0 200 400 600 


800 Vem? 1000 


Fig. 1. — Distribution in the atmosphere of fast nu- 
cleons (E > 3 GeV). 


B = 0.89. 


3:2. — Numerical calculations. 


Fig. 1 gives the absor- 
ption .curves of the total 
and of 
its protonic and neutronic 
components. The absolute 
intensities depend on the 


nucleonic radiation 


(24) This assumption is not verified experimentally for artificial mesons of low 
energy. However the theoretical arguments which give a greater probability for the 
production of non-identical particles do not appear to hold for the very fast nucleons 


considered here. 


(25) I. F. Quercıa, B. Rispoti and S. Scruri: Nuovo Cimento, 7, 715 (1950). 
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choice of values for H, and P, discussed in 45. The distribution arising 
from primary protons is represented by the dotted curve. The maximum 
which occurs there is absent from the curve for all neutrons owing to the 
large number of them contained in «-particles and heavier nuclei of the primary 
radiation. The slight maximum that the experimental curve seems to show 
at about 20 g/cm? is not a difficulty for our theory since we have completely 
neglected slow neutrons. | 


4, — Analisis of the mesonie component. 


4.1. — Creation of mesons. 


The scheme of n. 2 easily yields the espressions for the numbers dy" (W;,1l)4W;, 
dıg (W-,)AW., dig(W;,l)aW, of mesons rn, 7 and r° with total energy 
between W,, W. + dW, created between 1 ad 1 + di: 


dgt dW. = E mil. ml p(B, l)+ î °n(E, a} dE : dl, 


J 
\ 
dosdWe— | É + (1 — mel n(#, 1) + = (pE, | dE : di, 
J 
| 


1 | 
dy AW, = : 20 — ml p(E, 1) + lew fic nl n(E, nag te 


Introducing the expressions (8) for p(B, 1) and n(E,1) one obtains 


dW, dl 
et ee 


& (4) va, exp (— xl) + 3P,(— Pl)] 


1 JAN aW.di 
(9) aug AW = 5 # (Te [7H, exp (— al) — SP BI) 


£* 


5 We Nasa dw,, dl 
dig° AW; = 12 sH, (=) exp (— al) ve 


We may assume that all the rn’, 7 meson quickly (~ 107$ s) (**) transfornr 
into m+ mesons according to the scheme 


mt > rt + n° (n°, probably a neutrino), 


. (2) The processes of interaction of x-mesons with nuclei is negligible in compa- 
rison with that of disintegration in non-condensed matter. 
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since we do not consider particles with very high energy. The numbers 
dg=dW., give thus the numbers of tx mesons created in the length di. By 
putting q = 0.8e* © 0.2, we get: 


24° 


1 (we 1 
(9') dig=E(W, )AW= (=) (7H, exp (— al) + 3P, exp (— fl) u > 
Recent experiments (27) on artificial and cosmic rays mesons prove that 
the neutretto r° almost instantaneously (the mean life is smaller than 107% 8) 
disintegrates into two photons following the scheme 


To > hy, + ho, , 


suggested (28) by FINKELSTEIN and by VERDE. Therefore we assume that 
neutrettos do not contribute to the meson component, but only to the elec- 
tronic one (?°). 


4. — Distribution of mesons in the atmoshpere. 
u+-mesons can disintegrate (mean life t 2 2.15-10-*s) into one electron 


(+ or —) and into two neutral particles (apparently neutrinos): this would 
decrease the number of + mesons in a lenght di by: 


ig ae 
(10) d; g=(W, l). dW = g=(W, I) —dW= 
T 
_ A (W, D dtdW= À g#(W, ) dW dl 
— mW q , At —— WI g ’ 2) 
where 
: 1 : 2 
Ba Inc 21 3 Gever 
ToCy 


Furthermore the ionization losses are given by 


J ga W W = —— zer 


| dW dl, 
df 


(7) R. E. MARSHAK: Phys. Rev., 76, 1736 (1949); A. Lovati, A. Mura, G. SALVINI 
and G. TAGLIAFERRI: Nuovo Cimento, 7, 943 (1950); M. F. Kapton, H. L. BRADT and 
B. PETERS: Helv. Phys. Acta, 28, 28 (1950): R. BJORKLUND, W. E. CRANDALL, B. J. 
Moyer and H. F. York: Phys. Rev., 77, 213 (1950); W. K. H. PANOrSKY, LEE AAMODT 
and H. F. York: Phys. Rev., 78, 825 (1950); A. G. CARLSON, J. E. Hooper and D. T. 
Kine: Phyl. Mag., 41, 701 (1950). 

CS) R. FINKELSTEIN: Phys. Rev., 72, 415 (1947); M. VERDE: Ric. Scient., 17, 
2031 (1947). 

(2) H. W, Lewis, J. R. OPPENHEIMER and 8. A. WOUTHUISEN: Phys. Rev., 73, 
127 (1948). : 
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where we assume (°°) 


a 
dl 


| — 6 — const 2 1.48 : 108 eV. 
J 


The balance will yield 


DJS 


12 — 
(12) ol 24 q 


W\-s-1 D Se d Se 
ai [7H, exp (— al) + 3P, exp (-Bl)]— B = 48 a 
The solution is 


(13) g=(W, 1) = RW, 1) f Be BW, D) + SDW, n 


where 


? 


psy B|(W + 5) 
(14) Bm (2) 122) 


dx : 


(15) BW, 1) = | exp (— re +00 (3 HI 


age +6) —s—1 
0 


These results hold for W > W,—qH.. If W< Wo, equation (12) reads 


a Vays py 


0 Dre Wi ras 


since no mesons of such energy are created directly. The solutions of (16) 
which join smoothly the solution (13) (valid for W> W,) are 


(17) gt(W, D) = R(W, 1) |7 = @, (Wo, 1) + 3®p(Wo, |. 


The expressions for the total mesonic component are 


H H 
g(W, D = 14 R(W, DBA(W,D , g<(W, 1) =145 RW, DBa(Wo; D : 
0 0 
These results hold for magnetic latitudes higher than the latitude y for which 
the cut-off energy for protons H, equals the critical energy E, for the pro- 


duction of mesons (1°2). 


() The ionization loss is thought to be about 2.08: 107 eV for one cm of Hg. 
Considering o = 2.27: 10% em?, 1 = 1 corresponds to about 72 g/em?, thus giving 
the assumed value of à. 


2 - IL Nuovo Cimento. 


18 P. CALDIROLA, R. FIESCHI and P. GULMANELLI 


43. — Positive excess of the mesonie component. 


The positive excess is à good check on any theory on the mechanism of 
meson production. The usual measure of the positive excess is 


gi(W, ) — g-(W, 2) 
ME = 
area 


the differential positive excess so defined is a function of the energy of mesons 
and of the altitude. Our theory yields, by using (13) (16), 


for W>W, MM 


and for W<W, e(W, 1) = - = 


The positive excess for W< W, should be constant: obviously this result 
derives from assuming o constant for H > E, and zero for E< E,. Actually o 
will certainly not vanish for E< E,, but will go to zero gradually. So our 
value of e(W,l) is an upper bound for W< W,: a lower bound would be 
obtained by applying in this energy interval the formulas valid for W > W,. 


44. — Effects of the multiple production of mesons. 


In the collision of nucleons of high energy (above some 10 GeV) production 
of mesons with multiplieity higher than one can occur. These processes are 
less important than those with multiplicity one for cosmic rays in the atmo- 
sphere: however they might exert an appreciable influence on some quantities 
(e.g. the positive excess of the mesonic component). 

To evaluate the importance of multiple processes, we have computed the 
influence on the distribution of mesons of the following production with 


multiplicity 2: 


| P+ P+7°+ 7° | N+N+7°+ 7° 
ai tt - Nee + = 
P4p= P+ Pto3mr+_14 ru N+rt+r 
P+ N+7°4 nt N+P4+ r°4+ 7 
N+N+nt-+ nt | P+P+r +r 


[ P+N+m+n0 

TR i a, 

NEP |) Ppp yew 
| N+N+ 79+ xt 
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Equal probability is assigned to the various possibilities in each collision process. 
One obtains the following equations for the proton and neutron distributions 
in the atmosphere: 


dP(E, |) 9 E af E 
DR RE SSD Dee qe À It] roses phen) WE ! Dee) t 
PE DS] PUD pad à] NW, DAE 
(18) ro CE 
DN(E, |) 9 [ E to ff E 
Be NE) he ae D INA 
- EDS | ND zz, | PUY Ge a 
Ele Ble 
their solutions are: 
1 268 
P(E,) =; E-‘(H, exp (— al) + P, exp (— Bl)] | a jet — 0.574 
2 8 + 1 
(19) : where 
N(E, 1) = 2 E-[H, exp (— al) — Py exp (— fl) | Le =1 en — 0.947: 


The corresponding equations for the distribution of mesons produced in double 
processes are: 


[ 25 (W, 1) ak (7) (9%. exp (— al) + 3P, exp (— fl) |— 


Me as 


qd | 
— pF (WD | og (W, I) PAU 
(20) B Wi + Ô Ww ; for We 2 , 
dg (W, 1) gz (W, 1) ogi (W, 1) Wo 
L 2 — a Ka à A 
N a a 
Their solutions read: 
3 H W 
| g# (Wl) = 55 RW, D) Ê ms D(W, 1) + 384 (W, 0) , W> oe 
4 0 
W 


| 0=(W,) = — R(W, 1) Ê 7 D. (5 ) + 30», (= ; ) ; Wen ; 
It is easy to verify that the positive excess ealeulated on the ground of meson 
production of multiplieity 2 is much smaller than the one for multiplieity 1. 
On the other hand the probability of the «double » processes increases with 
energy. We shall assume that in a certain interval of energy no processes 
of multiplieity greater than 2 occur: then the number g(W,1) of mesons of 
energy W observed at a depth 7 can be written as a sum of the number of me- 
sons created, with the probability a(W,1), in double processes and of the 
number of mesons created, with the probability 1 — a(W, 1), in single processes 


g(W, D = gt(W, 0) + g(W,) = 1 — aW, Dlg (WY) + a(W, 1)g.(W, 1). 
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The differential positive excess is given by 


> g'(W,D— g-(W, ı) 
21 (W,)=2% a Se 
or 
_ 5 ba, Dilgt (W, D — 9 (W, D+ aW, Digs (W, 0 — gs (W, D] 
[1 — a(W, )]g(W, D + a(W, DgAW, 1) ; 


where the function a(W, 1) 
can be determined from 
experimental results (see 
the following number). 


45, — Numerical caleulations 
and comparison with 
experiments. 


Fig. 2 gives the curves 
for spectral distribution of 
the mesonic component for 
1 = 14.4 (sea level) and for 
l= 5 (about 8 000 m a. s.1.): 
H. has been taken equal 
to 3 GeV to have a latitude 
effect beginning at 9 = 45° 
in agreement with experi- 
ment. The curves are very 


PRE En hae ET UC or 3, . Similar to the experimental 
P(MeV£) ones obtainable by extra- 

Fig. 2. — Spectral distribution of mesons at different polating the results of 
dephts/(at sea level and at 360 g/ em?). HALL (#1) (experiments per- 

© Wilson formed at Mount Evans, 


about 4350 m) and of 
SCHEIN and coworkers (3?) (experiments performed at 5-6000 m) for the high 
altitude part and using the results of WILSON (3?) and BLACKETT (*) for the 
low altitude part. The comparison of the two curves shows also the well 
established increase of the number of slow mesons with altitude. 


(51) D. B. HALL: Phys. Rev., 66, 321 (1944). 
(2) M. ScHEIN, O. E. WoLLAN and G. GROTZINGER: Phys. Rev., 58, 1027 (1940). 
(33) J. G. Wırson: Nature, 158, 414 (1946). 


G 
(34) P. M. S. BLACKETT: Proc. Roy. Soc., A159, 1 (1937); see L. LÆPRINCE-RIN- 
GUET and J. CRUSSARD: Journ. de Phys. et Rad., 8, 207 (1937). 
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Numerical integration on the energy spectrum of mesons for different values 
of 1 has yielded the number of mesons at different altitudes. Mesons have 
been classified as fast (W > 0.33 GeV) and as slow (W< 0.33 GeV) to com- 
pare our results with the experimental curves of Rossi (%). Fig.s 3 and 4 
show the results obtained. H,, the only parameter yet indetermined, has 
been chosen to give the cor- 
rect number of fast mesons "FT | | | De I ware 
at sea level: the data on the ® 
distribution of fast mesons are 3 
more exact than those for  § 

S 
> 


a 


a 


slow mesons. Fig. 3 shows 
that the theoretical’ absorpt- 
ion is slightly steeper than the 
experimental one: the posi- 10° 
tion of the maximum is at 

200 g/em?, in good agreement 5 F 
with experiments. The agree- 
ment is certainly satisfact- 


JE 


ory, if one considers the sim- 2 


plifying assumptions of the 


theory (one-dimensional ap- eh 
. . . 4 0) 200 400 600 800 9m? 1000 
proximation, assumption on cm 


the cross section, etc.) and Fig. 3. — Distribution of fast mesons in the 
atmosphere. 


the uncertainty of the expe- 
rimental data. 

The value found for A, 
(2.8-101% em~*s~'sterad eV.) 
has been used to ealeulate the distribution of slow mesons: fig. 4 gives the re- 
sults. The number of slow mesons is slightly smaller than the experimental 
one: however the simplifying assumptions on the’ cross Section reduce appre- 
ciably the probability of produetion of slow mesons and, on the other side, 
experimental data are quite uncertain. 

The results obtained allow us to determine the parameter A which measures 
efficiency of production of mesons by the nucleons contained in «-partieles 
and heavier nuclei. Indeed the experimental results (3°) on the spectrum 
and the nature of particles (protons and nuclei) arriving at the top of the 


—___— Experimental curve (Rossi) 
— — — — Theoretical curve. 


(35) B. Rossi: Rev. Mod. Phys., 20, 537 (1948). 

(36) The observations on the particles incident at. the top of the atmosphere can 
be interpreted by assuming an intégral spectrum of the. type 2.4: 1016. E28, and 
considering the protons to be a bit more than 70% of the incident particles (equal to 


about 30%, of the total number of nucleons). 
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atmosphere give the values 


P, = 1.68 - 101% em? gs“! sterad—1 eV!-® 
H, = 2.84 : 1016 em? 87? sterad-! eV-18: 


since equation (4) holds, the value already found for H, gives A — 0.39. 

This result is very useful 
in interpreting the experi- 
mental data on the positive 
excess of the mesonic com- 
ponent. The existence of this 
excess was proved experimen- 
tally by various authors (°”) 
and particularly by the direct 
experiments of BERNARDINI 
and coworkers (35): these mea- 
surements seem to ‚indicate 
a value of 15 (+5) % for 
the integral excess (cut-off 
energy 0.3 — 0.5 GeV) at sea 
level. The first measure- 


ment (*) of the positive ex- 


cess of mesons for given e- 
0 200 400 600 800 G/.m2 tooo Nergy showed à rather un- 
certain decrease of the diffe- 
rential excess with increasing 
—— Experimental curve (Rossi) energy: later results of Bro- 

— —— — Theoretical curve. DE (*), of the group of Pa- 

dova (41), of WILSON and 

OWEN (42) and of CARO, PARRY and RATHGEBER (*) show a noticeable in- 
crease of the positive excess with energy up to 5 GeV while, at higher 


Fig. 4. — Distribution of slow mesons in the 
atmosphere. 


(37) H. Jones: Rev. Mod. Phys., 11, 235 (1939); D. J. HuGKes: Phys. Rev., 57, 
592 (1940). 

(38) G. BERNARDINI, M. Conversi, E. PANCINI, E. Scrocco and G. C. Wick: 
Phys. Rev., 68, 109 (1945). 

(%) M. Conversi, E. PAncını and O. Pıccioni: Nuovo Cimento, 3, 372 (1946); 
N. NERESON: Phys. Rev., 73, 365 (1948). 

(4) R. B. BRODE: Bull. Am. Phys. Soc., 5, 24 (1949); Como Confer. Sept. 1949. 

(#1) P. Bassi, E. CLEMENTEL, I. FrLosoro and G. Puppr: Como Confer. Sept. 1949. 

(2) J. G. Wırsox: Como Confer. Sept. 1949; B. G. OwEn and J. G. WıLson: Proc. 
Phys. Soc.. A 64, 417 (1951). 

(4) D. E. Caro, J. K. Parry, and A. D. RATHGEBER: Nature, 165, 689 (1950). 
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energies, the accurate measurements of Wrrsox and OWEN show a decrease 
of the excess. 

Fig. 5 shows the theoretical curves for the positive excess due to single 
and double meson production and also the experimental curves of WILSON 
and Owen. Below 5 GeV the first theoretical curve agrees well with the 


025 


020 


015 


010 


0 2 4 6 8 10 p(GeV/) 
Di, Hg = Differential positive excess at sea level. 
A Wilson. 


experimental one, indieating that mesons in this energy interval have been 
created by nucleons (with an energy <25 GeV) by single processes. The 
decrease of the positive excess of mesons of energy higher than 5 GeV indi- 
cates instead that nucleons of energy higher that 25 GeV create mesons by 
multiple processes with multiplicity increasing with energy. Assuming that 
only single and double processes occur Up to a certain energy, formula (21) 
holds. Up to energies of 11 GeV, it is possible to determine the function a(W) 
in it to obtain a fair agreement between theoretical and experimental results; 
above 11 GeV this is not possible any more, indicating that at these energies 
multiplicities higher than two are present. Since the average energy W of 
a meson is connected with the energy of a generating primary by the relation 


W 
pe tae) RU os an 
q q 


we can obtain the dependence of the average multiplicity on the energy of 
the incident nucleon 


K(E) = 2a + (1—a)=a+i. 
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This is illustrated in fig. 6 (*). The calculations seem to show that in the col- 
lision of two nucleons the production of mesons is practically single up to 

energies of 25 GeV; while, 
2 Er for energies between 25 

and 70 GeV, double pro- 
duction is important and 
is actually prevalent be- 


15 yond a certain point. For 
energies above — 70 GeV 
the multiplieity is higher. 
Fig. 7 shows that on 
E(GeV) passing from I = 14.4 (sea 
= 5 (360 g/cm?) 
Fig. 6. — Dependence of the production multiplicity on the positive excess does 
the energy of the incident nucleon. slightly increase. This re- 


20 40 60 80 100 level) ; 


sult of our calculation 
agrees even quantitatively with the measurements of QUERCIA, RISPOLI and 
SCIUTI (4°) at high altitudes. Even more significant is the agreement of our 
theory with the results of these authors (#) on the dependence of the po- 
sitive excess on the direction of the mesons. 


915 
0 2 4 6 8 P(GeVÆ) 
Fig. 7. — Differential positive excess at sea level and at 360 g/em2. 


(*4) In determining the function K(E) we have not considered the ionization losses 
of mesons, from their production quota to sea level. Owing to this simplification phe 
values of F, in the figure corresponding to K, are to low of a few GeV. 

(5) I. F. QuercrA, B. RıspoLı and S. Scıurı: Nuovo Cimento: 5, 397 (1948); see 
M. Conversi: Phys. Rev., 76, 311 (1949). 

(78) I. E. QUERCIA, B. RıspoLı and 8. Scıutı: Nuovo Cimento, 7. 715 (1950). 
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All our calculations have been done on the assumption of a very parti- 
cular mechanism of energy transfer from the incident nucleon to the scattered 
one and to the mesons produced. Though our results depend on the assump- 
tions made, we believe that reasonable modification of them will not produce 
substantial changes. Furthermore, since experimental data are uncertain, our 
considerations are purely indicative. 


5. — The electronic component. 


5'1. — General remarks. 


The electromagnetic component of cosmic rays is generated by different 
processes: the most important is now considered to be the disintegration of 
the nentretto into two photons. 

Since it appears plausible that the number of neutrettos produced in the 
interaction between nucleons is comparable with the number of charged me- 
sons, it is obvious that this process should be the main source of the electro- 
magnetic component, expecially at high altitudes. Important contributing 
effects, however, are the disintegration of „-mesons and their knock-on col- 
lisions on the electrons of the air. On the other hand, the knock-on of protons, 
the direct production of y-rays in stars or in collisions between nucleons and 
the absorption of x mesons by nucleons, according (*) to BRUNO and to 
MARSHAK, do not give an appreciable contribution and are not cousidered in 
our calculations. 


52, — Disintegration of the neutretto. 


The assumption of isotropie emission of the two photons in the center of 
mass system gives the value W,/2 for the average energy of each photon. 
The number of photons produced between I and 1 + di is then given by 


dıy(W, l)dW= 4d,9°(2 W, I) OU 


Assuming that 7° disintegrate immediately after creation, equation (9) 
gives the following average spectrum of formation: 
N Aa nl 
dl, DAW = =~ Sholay exp (— al) 
= 


dw 
D = 


=O). 
x: 


(47) R. E. MARSHAK and A. S. WIGHTMAN: Phys. Rev., 76. 114 (1949); B. BRUNO: 
Ark. Phys., 1, 19 (1949). 
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Naturally this holds only for latitudes such that E,„< E,; then the mi- 
nimum value of W in the spectrum of photons is e*Æ,/2. 

Under these boundary conditions the theory of showers (4%) in Arley’s (*) 
approximation yields the following electronic integral-spectrum 


[ Bo |B; exp! Ai(s)t}—-exp[—aad] exp [A,(s)t]—exp Fe) 
HUE PRE A,(s) ao A,(8)-- aw Fe 
29 exp [A,(s)t] 1— exp [— À(s)o] exp [- aat] | 1 — exp [axo] 
122) ae Aiea A,(s) + ax ac A,(8) + ax 
| exp [A (s)t] .1—exp[— À(s)o] | exp [— aut] 1— ul 
(8) B)+ax | ax As(s) + ax JJ 


For details see paper IV. 


5°3. — Electrons from w-mesons disintegration and knock-on electrons. 


It is almost certain (5) that w= mesons disintegrate in one electron 
(+ and —) and two neutrinos. The spectra of the decay electrons give for 
the mean energy E’ in the reference system of the meson a value between 
35 and 40 MeV: the most probable value is 36 MeV. In the laboratory system 
E = 36 ( W/u) — 4 W; then the spectrum of formation of the decay elec- 
tron is: 


Bye, 


The problem of the shower originated from this source has been solved in 
paper IV by means of Mellin transforms and the saddle point method, both 
for energies higher and lower then the critical energy E, of the air. 

Similarly one could evaluate the contribution to the electromagnetic com- 
ponent from the knock-on of y-mesons: however the formula for the highest 
energy given to the electron shows that the contribution of slow mesons can 
be forgotten. 

The calculations for fast mesons do not need to be carried out, since the 


(4) B. Rossı and K. GREISEN: Rev. Mod. Phys., 18, 240 (1941). 

(#8) N. ArLEY: Proc. Roy. Soc., A 168, 519 (1938). 

(50) E. P. Hıncks and B. PontEcorvo: Phys. Rev., 75, 698 (1949); R. B. LEIGHTON, 
€. D. ANDERSON and A. J. SERIFF? Phys. Rev., 75, 1432 (1949); J. D. Davızs, W. O. 
Lock and H. Murrunap: Phyl. Mag., 40, 1250 (1950). 
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ratio of the intensity of the 
electromagnetic component 


for knock-on to those for de- 
cay, for the same spectrum of 
generating mesons, is already 
known (51). 


54. — Numerical valuations. 

Fig. 8 compares our re- 
sults with the measurements 
of POMERANTZ (5?) and with 
the curve of Rossı (#): the 
agreement between theory 
and experiment is satisfact- 
ory. The dotted curve gives 
the partial contribution of 
neutrettos: the depth corre- 


sponding to the maximum 
of the curve agrees with ex- 


periment and its position is 
not very sensitive to the va- 
lues chosen for the parame- 
ters; the intensity of the curve 
near the maximum almost e- 
quals the total intensity and 


the difference of the two lo- 


gically increases towards low 
altitude. 
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Fig. 8. — Distribution of the electronic component 
in the atmosphere. 
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The calculations show that slow mesons are much less important than 


fast ones. 


(1) B. Rossı and S. J. KLAPMAN: 
B. N. CaccrapuorTi and B. QUERZOLI: 


Phys. Rev., 61, 414 (1942); G. BERNARDINI, 
Nuovo Cimento, 3, 349 (1946); L. JANossy: 


Cosmic Rays (Oxford, 1948), page 247. | 
(5) M. A. Pomprantz: Phys. Rev. 75. 69 and 1721 (1949); 77, 830 (1950). 


28 P. CALDIROLA, R. FIESCHI and P. GULMANELLI 


6. — Latitude effects. 


6'1. — Latitude effect of the nucleonie component. 


We have to integrate equations (6’) of the nucleonic shower under proper 
initial conditions. Let H, > H, be the cut-off energy of the earth’s field at 
the magnetic latitude m: the initial conditions read 


* 
0 


ER 

P(E, 0) = (Po + EX EB (W<E<E,) 
H* | 

P(E, 0) = (r, ET 

NEO 


Jake 
NEO) = AB a (y <£). 


——,— 
* 


A detailed study of equations (6’) under these conditions was carried out im 
paper II. ÆA(E,1) = P(E,1) + N(E,D and K(E,1) = P(E,1)— N(E, 1) sa- 


tisfy the equations: 


[ H(E / 
el en Sen | H(E',) 
ol € x IE 
(23) | Eie 
A E b n LA 
le = Dae a, | en ee 
Ele 


with the new boundary conditions 


j HMO) HE) (O< B= Ke) [K(#, 0) = PES Ot ae) 


(24) 
| H(#,0)=H,E-- (E,<E); | K(E,0) = PE" (E,<E). 


Therefore we simply have to study equations of the type: 


F(x, y) on |. d 
25 es , = M(E 
(25) x + F(x, y) = | P(E, y) £ 


with boundary conditions 


12) AA (<a) IC) Zaren). 
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The solution F(x,1) can be written in polynomial form 


(2 > Xp) 


Bla) Am exp (, ae = 1) Y 
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Die © ET = Te + exp (arm) | Er > ee 
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(er < 2 < eo; v 21). 
0 5 10 Ey (GeV) 15 


The results of the numerical calculations 4 
are given in figs. 9, 10, 11. Fig. 9 re- 
presents the ratios 


080 


| ner 
Hy(Bx; t 070 
AIRE fe = a) 


Be) 

between the number of fast nucleons at a 
the latitude g and the number of nu- 
cleons of the same type at latitudes 
greater than % (45°) where the effect  oso 
of the cut-off is not felt any more. The Fig. 9. — Ratios between the total 
curves have been calculated for I, = 14.4 number of nucleons with KH > EB, at 
and l,=5. At sea level the ratios are latitude p and at latitudes > 45°. 
practically the same for the various types 
of nucleons, but this is not so at higher altitudes where the values of Ry be- 
come appreciably greater than those of Rp. 

Fig. 10 represents, for I, = 14.4, the ratios 
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between the distribution functions P,(E, lo), Ny(E, lo), H,(£, lo) for the values 
5, 10 and 15 GeV cf the cut-off energy and the corresponding distributions 
without cut-off. Fig. 11 represents the same ratios for 1 = 5. 


E(GeV) ıs 
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096 


Fig. 10. — Ratio between the integral spectrum of nucleons at latitude gy and the 
corresponding spectrum at latitudes > 45° (1, = 14.4). 
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Fig. 11. — Ratio between the integral spectrum of nucleons at latitude g and the 
corresponding spectrum at latitudes > 45° (1, = 5). 
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The curves of fig.s 10 and 11 have been used to study the latitude effect 
of the mesonic component. The magnitude of the effect at sea level is 4% 
at 0°. The effect increases appreciably with the altitude and at ~ 8000 m 
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is 40% for the neutrons and 50% for the protons: this implies an increase of 
the mean free path for absorption, in passing from the northern latitudes to 
the equator, greater for neutrons than for protons. Fig. 1 shows that the 
absorption eurve for fast neutrons differs from the exponential, at depths. 
smaller than 300 g/em?, when no cut-off is present: at the equator this will 
oceur even at greater dephts. 

These results are in qualitative agreement with experiment (*), but no 
careful check is possible since the exponential data are as yet few and rather- 
uncertain: the available data on slow neutrons cannot be used in this con- 
nection. 

APPAPILLAI and MAILVAGANAM (5) did not find at sea level and g = 3° 
any latitude effect on the generating component of hard showers; but the 
uncertainty of their data can mask an effect of 4%. Rossi et al. (55) have 
found at 20° and 300 g/em? an effect of ~ 20% for protons of energy > 5 GeV 
and of ~50% for protons of energy > 0.4 GeV. 

We should conclude that the latitude effect that we predict for the nucleonic 
component is rather high: the same holds if we compare our results with the 
experimental findings of WINCKLER and coworkers (5%) and of WHYTE (5). 
Other theories and particularly MESSEL’s one give the same result: MESSEL ($) 
suggested that the exponent chosen for the primary spectrum might be too 
high. Recent experiments confirm this. 


62. — Latitude effect of the mesonie component. 


The latitude dependence of the spectral and of the absorption curves and 
of the positive excess of mesons is a good check on any scheme for the pro- 
duction of mesons. The utility of this check has been proved by various 
authors and particularly by M. FERRETTI SFORZINI (5%) who has shown that 
the latitude effect of the meson component is incompatible with high multi- 
plicity productions. 

To carry out the calculations we start from the formulas for the distri- 
bution in the atmosphere of the various nucleonie components (protonic, 
neutronic and total) at different latitudes (see preceding number) and from 


(83) J. STAKER: Phys. Rev., 80, 52 (1950); J. A. SIMPSON and E. HUNGERFORD: 
Phys. Rev., 17, 847 (1950). 

(54) V. APPAPILLAI and A. W. Mairvacanam: Nature, 162, 887 (1948). 
(55) A. J. Mc Manon. B. Rossı and W. F. Burpirr: Phys. Rev., 80, 157 (1950). 
(56) J. R. WINKLER, T. STIK K. Dwicnt, and R. Sasın: 79, 656 (1950). 
(57) G. N. Wnyre: Phys. Rev., 78, 321 (1950). 
(58) M. FERRETTI SFORZINI: Nuovo Oimento. 7, 196 (1950), 8, 909 (1951). 
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the general formulas for the distribution of mesons given in 4 One obtains: 


W\2/01+W) N R L ; a n = 
g=(W, 2) = | | eBil+W)W + öll— a) PGI). Sfr, W+ Al 


æ)| dx, 


where : 


OST EE i (> eee ~= (5: Say 
S- amd: Choy = 12? | 5") À dl, Gl = 6 a n ji dl. 


: Two arbitrary constants are contained in these formulas: the fraction e* 
of the energy of the incident nucleon transferred to the mesons created in the 
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Fig. 12. — Differential mesonic spectrum at sea level at different latitudes. 


collision with another nucleon; and the critical energy ÆZ, for the production 
of mesons. The calculations which do not involve the geomagnetic effects 
are fairly insensitive to the numerical values of these parameters; instead 
the results for the geomagnetic effect for mesons show a critical dependence 
on the values of «* and Æ, and allow one to determine them rather accu- 
rately. The best results are obtained assuming H, — 3 GeV, &* — 0.25 and 
keeping 125 g/cm? as the value of the absorption coefficient for the nucleonic 
component. x 

Fig. 12 gives the computed energy spectra of mesons at sea level and in 
the vertical direetion corresponding to latitudes identified by different values 
of the cut-off energy H,: the latitude effect appears most relevant for mesons 
of energy smaller than a few GeV. 
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Fig. 13. — Latitude effect of the mesonie component (at sea level). 


Fig. 13 gives the percentual latitude effect on the total number of mesons 
(energy > 0.2 GeV) obtained from the computed spectra at the various lati- 
tudes using VALLARTA’s theory (5) for the correspondence between cut-off 
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Fig. 14. — Differential positive excess at different latitudes (at sea level). 


(5) D. J. X. Montcompry: Cosmic Ray Physics (Princeton, 1949), page 55. 
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energy and magnetic latitudes. The theoretical curve has the characteristic 
behaviour of the experimental curve and the numerical values agree farly 
well (5). 

Fig. 14 gives the curves for the differential positive excess of the mesonic 
component at sea level for different values of the cut-off energy: the results 
are normalized to keep account of x-particles and heavier nuclei in the pri- 
mary radiation (see 45). The magnetic field lowers the positive excess for 
low energy mesons. 

Most experimental data as yet available give the latitude effect on the 
integral intensity of the mesonic component and not the dependence of the 
effect on the energy of mesons. A strict comparison of energy spectra of 
mesons at different latitudes determined under the same experimental con- 
ditions, as well as the comparison of the curves for the differential positive 
excess at different latitudes, could be of great use in checking the adequacy 
of our scheme as well as of any theory on the mechanism of production of the 
mesonie component. Also the study at different latitudes of the dependence 
of the energy spectrum and of the curve of the differential positive excess on 
meson direction would yield interesting data. 


Conclusions. 


We have seen that the measurements on the distribution of the different 
components of cosmic rays in the atmosphere are compatible with what one 
knows on the composition of the primary radiation and on the processes 
through which the nucleonic component transforms into the mesonie and 
electromagnetic components. 

The good agreement with experiment of the curves that we have computed 
with simple assumptions shows that the current ideas on cosmic ray processes. 
are substantially correct: the phenomena not considered (e.g. V particles pro- 
duction) and the eventual processes as yet unknown are not very important 
in the total energy balance. Our results allow us also to draw some con- 
clusion on the interaction between particles; in particular our scheme agrees 
well with the hypothesis of single meson production up to energies of a few 
times 10 GeV and of multiplicity slowly increasing with increasing energy of 
the interacting nucleons, in fairly good agreement with FERMI’s theory. 


(°°) D. J. X. MONTGOMERY: Cosmic Ray Physics (Princeton, 1949), page 127. 
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RIASSUNTO 


Partendo dai risultati sperimentali relativi alla composizione della radiazione cosmica 
primaria, viene sviluppata una teoria fenomenologica di tutti i processi che hanno luogo 
attraverso l’atmosfera. In particolare si calcola: a) la distribuzione della compo- 
nente nucleonica (protonica e neutronica) veloce; b) la distribuzione della componente 
mesonica, il suo spettro energetico e il suo eccesso positivo a diverse altezze; e) la distri- 
buzione della componente elettronica. Vengono poi studiati accuratamente gli effetti 
di latitudine dovuti al taglio del campo magnetico terrestre. Tutti i risultati teorici 
ottenuti sono confrontati con i più recenti dati sperimentali. 


IL NUOVO CIMENTO VOL IX EN 1 Gennaio 1952 


Luminescenza di soluzioni per effetto di raggi y 
rivelata mediante un fotomoltiplicatore 931 A. 


A. CICCONE 


Istituto di Fisica dell Universita - Pisa 


(ricevuto il 25 Ottobre 1951) 


Riassunto. — Nella presente nota si riferisce su esperienze eseguite al 
fine di studiare la luminescenza di soluzioni per effetto di raggi y, e 
Vintensita della luminescenza al variare della concentrazione; si ripor- 
tano due possibili spiegazioni qualitative del fenomeno della luminescenza. 


1. — Introduzione. 


Molti liquidi organici, quando sono allo stato di grande purezza, eccitati 
con radiazioni di alta energia (raggi «, 8, y, neutroni, ecc.) emettono una debo- 
lissima luminescenza che diventa assai più sensibile se nel liquido sono disciolte 
piccolissime dosi di speciali sostanze come antracene, fenantrene, naftalina, 
terfenile, carbazolo, ece. 

Questa luminescenza dipende dalla specie della sostanza in soluzione, dalla 
sua concentrazione e dalla natura del solvente adoperato. Per lo studio di tali 
luminescenze male si prestano i soliti metodi spettroscopici o spettrografiei a 
causa della debolissima intensità di esse. Ottimamente invece si prestano i 
fotomoltiplicatori per la loro grandissima sensibilita e per il carattere quanti- 
tativo dei risultati che se ne possono ricavare. 

AGENO e collaboratori (1) hanno studiato la luminescenza prodotta dai 
raggi y in una soluzione di naftalina in xilolo sia in funzione della concentra- 


1 
zione, sia in funzione dello spessore della soluzione; REYNOLD (7), KALLMANN (?) 


() M. AGENO, M. Cutozzorro e R. QUERZOLI: Acc. dei Lincei, 6, 626 (1949). 
(2) G. REYNOLD: Nucleonics, 6 ,69 (1950). 
(?) H. KALLMANN: Phys. Rev., 78, 621 (1950). 
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hanno esaminato la luminescenza di liquidi organici puri, e di soluzioni ecci- 
tati con diverse sorgenti di alta energia; REYNOLD, HARRISON, SALVINI (%) 
hanno studiato la luminescenza di liquidi organici puri, di soluzioni in deter- 
minate concentrazioni e di liquidi organiei in cui erano state disciolte quan- 
tità diverse di due sostanze luminescenti. I numeri che questi autori ripor- 
tano, rappresentano cid che essi definiscono « efficiency » di questi liquidi. 

Pit organico & il lavoro di KALLMANN e FURST (5) i quali hanno studiato 
sia la luminescenza di soluzioni diluitissime di antracene, carbazolo e terfenile 
in xilolo in funzione della concentrazione, sia i] comportamento deJl’intensita 
della luminescenza di aleune soluzioni per aggiunta di un’altra sostanza lumi- 
nescente. Questi autori suggeriscono inoltre un’interpretazione dei fenomeni 
osservati basata sul trasporto dell’energia di eccitazione dalle molecole del 
solvente alle molecole del soluto. 

E da notare, peraltro, che la maggior parte delle misure eseguite in questo 
campo (6) non sono facilmente confrontabili le une con le altre. 

In questo articolo si riportano i primi risultati di una serie di misure intese 
a investigare in modo sistematico e, per quanto possibile, esauriente le carat- 
teristiche di liquidi luminescenti. I risultati qui riportati si riferiscono a solu- 
zioni in xilolo di aleune sostanze che ho potuto procurarmi in ottime condi- 
zioni di purezza. Per questo studio mi è stato prezioso l’appoggio dell’Istituto 
di Chimica Farmaceutica della nostra Università (*). 


2. — Dispositivo sperimentale. 


Tra i vari tipi di fotomoltiplieatori i pit adatti per il nostro uso sono 
l’1P28, il 931A, 171 P22, i quali differiscono fra di loro perchè utilizzano super- 
ficie fotocatodiche ent curve di risposta spettroscopica diverse l’una dal- 
l’altra; io ho potuto avere a mia disposizione un 931A. Questo ha una super- 
ficie fotocatodica di antimonio-cesio (?) e una zona di sensibilita spetbroscopica 
compresa fra 4000 e 5000 À con un massimo intorno à 4500 À. Esso servirà 
dunque, per rivelare luminescenze le cui Junghezze d’onda sono comprese in 
questo intervallo. 

Il primo problema & stato quello di mettere il fotomoltiplicatore (FM) in 
condizioni di massimo rendimento. È noto che per usare un FM si possono 
seguire due metodi: o la registrazione e il conteggio dei singoli impulsi, o la 


) G. Reynorp, F. B. HARRISON € CG. Sazvinr: Phys. Rev., 78, 488 (1950). 
5) H. KALLMANN e M. Furst: Phys. Rev., 79, 858 (1950); Nucleonies, T, 69 (1950). 
) Cfr. per es., G. REYNOLD: loc. cit. (2); G. REYNOLD e collaboratori: loc. eit, (*), 
*) Al Direttore prof. R. DE FAZr, all’assistente prof. CARBONT e al dott. SEGNINI, 
esprimo i miei più sentiti ringraziamenti. 

(7) J. W. Corrman: The Scintillation Counter (Proc. of the I.R.E.). 
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misura della corrente media mediante un galvanometro. Ho usato questo 
secondo metodo adoperando un galvanometro Deprez d’Arsonval la cui sen- 
sibilità, ad un metro della 

Gi Bel Rowe Rae Ber Ren uke Rag (Re scala dallo specchio, era di 
1,3. 0 PAS per sini: 

La fig. 1 rappresenta lo 
schema elettrico adoperato ; 
le tensioni erano mantenute 
costanti con appropriati sta- 
bilizzatori. 
fi fils Per avere un buon fun- 


© zionamento del FM con- 
2 viene proteggerlo da varia- 
à Ds zioni di temperatura; per 
= AVES questo motivo esso & stato 

ern foooy +1 racchiuso in una scatola di 


lamierino d’ottone, a doppia 
parete, con il fotocatodo in 
eorrispondenza di una pic- 
cola finestra rettangolare (praticata in una delle pareti) alla quale si accosta 
un diaframma in ottone pure sfinestrato in corrispondenza (fig. 2). 

Per evitare radiazioni estranee 
provenienti dall’ambiente, la scatola 
e stata racchiusa in una cassa di 
legno (90 x35x29 cm?, spessore 
em 3) annerita internamente. 

Le connessioni fra il FM e gli 
alimentatori sono state effettuate 
in modo da garantire un perfetto 
isolamento e da prevenire qualsiasi 
infiltrazione di luce. 

Come & noto il FM presenta, 
anche nella perfetta oscurità, una 
debole corrente «corrente di zero », : Fig. 2. 
la quale, se non si prendono spe- 
ciali precauzioni pud superare in intensità la corrente segnale da misurare. 

Per questo si prendono in considerazione le differenze d’intensitä, ma se 
si ricerca grande sensibilitä, bisogna, naturalmente, che la corrente di zero sia 
ridotta a un valore molto basso. Cio si otterrebbe facilmente raffreddando 
il FM con una miscela di ghiaccio e sale, o con anidride carbonica e acetone, 
ma riconobbi che questi metodi, per le mie esperienze, portavano ad altri 
inconvenienti e quindi ho preferito operare alla temperatura ambiente. 
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Varie prove dimostrarono che la corrente di zero, in buona parte era dovuta 
a dispersione per imperfetto isolamento, ma con opportuni accorgimenti si 
© riusciti a ridurre la corrente di zero a 0,3 - 10 A raffreddando il FM alla 
temperatura di zero gradi centigradi, a 2,4: 10° A alla temperatura ambiente 
di 13°C, e a 43-10 A alla 
temperatura estiva di 27 °C; 
in ogni esperimento la tensione 
fra un dinodo e il successivo 
era di 100 V. 

Una valutazione delle buone 
condizioni sperimentali cosi rea- 
lizzate, è stata ottenuta nel 
modo seguente: jnviando sul 
fotocatodo una radiazione lu- 
minosa d’intensitä determinata, 


si leggeva una deviazione al 
galvanometro; la differenza tra Fig. 3. 

questa lettura e quella corri- 

spondente alla corrente di zero esprimeva l’effetto della radiazione. Il rapporto 
tra questa differenza e la corrente di zero si pud assumere come una misura 
della sensibilita del FM (§). 

Questa luce d’intensità determinata e fissa era fornita da una lampadina 
alimentata da una corrente d’intensità molto piccola (tale da rendere il fila- 
mento appena luminoso) che veniva mantenuta esattamente costante. 

La lampadina era posta ad una distanza di 15 em da un diaframma collo- 
cato immediatamente davanti alla finestra della scatola (fig. 2) e avente un 
foro del diametro di 3 mm: essa era sistemata in modo che la sua luce cadesse 
perpendicolarmente sulla superficie del fotocatodo. 

Nei grafici della fig. 3 è riportata la sensibilità del FM in funzione della 
tensione fra i dinodi estremi (punti 4 e B di fig. 1) per due diverse temperature. 

Essi mostrano che al crescere della tensione la sensibilità eresce lentamente 
e suggeriscono una tensione di 100, 110 V, fra un dinodo e il successivo, come 
la pil opportuna. 


3. — Risultati sperimentali. 


Si sono preparate due soluzioni in xilolo, una di naftalina e una di ace- 
naftene, entrambe prodotti purissimi rispettivamente delle Case Merck e 
Fraenkel, alla concentrazione di 0,18 g/em*. Lo xilolo € stato, per maggiore 
precauzione, privato da eventuali tracce di umidità e ridistillato in laboratorio. 


(8) M. AGENO, M. Cutozzorro e R. QUERZOLI: Nuovo Cimento, 6, 621 (1949). 
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Le vaschette contenenti le soluzioni sono state appositamente da me co- 
struite con vetri a pareti sottilissime (circa 3/10 di mm) ed hanno una larghezza 
di mm 19 ed una altezza di cm 5; esse potevano adattarsi vieinissime alla 

finestra della scatola di ottone 
À x10? che conteneva il FM. 
® La sorgente, un sale di 
Radio di eirea un millicurie, 
è stato sistemato in una sca- 
toletta di alluminio dello spes- 
sore di 8/10 di mm, onde as- 
sorbire i raggi ß. Essa irra- 
diava le sostanze in esame 
attraverso un foro di 3 mm 
praticato in un blocchetto di 
piombo protettivo. 

Si è dapprima studiata la 
luminescenza di un cristallo 


2 # 
© di antracene molto traspa- 
rente, di g 3,76, il quale ve- 
Tr ua] . . 
700 800 900 1000 v niva collocato con una faccia 


Fig. 4. - a) cristallo di antracene; b) soluzione di  immediatamente davanti alla 
acenaftene in xilolo; €) soluz. di naftalina in xilolo. finestra del FM, mentre sulle 

rimanenti facce veniva rive- 
stito con un sottile foglio di stagnola per raccogliere anche la luce emessa 
nelle altre direzioni. 

Come misura de) potere luminescente del cristallo si à assunta la diffe- 
renza tra la corrente anodica misurata con e senza la sorgente radioattiva. 
Il grafico a della fig. 4 & stato ottenuto riportande il potere luminescente in 
funzione della tensione applicata tra A e B (fig. 1). Gli altri due grafici be ¢ 
si riferiscono rispettivamente alle soluzioni (7,520 8, concentrazione 0,18 &/em?) 
di acenaftene e di naftalina e sono stati riportati per confronto con il grafico a. 
Uguale geometria in tutti e tre i casi. 

AJ fine di verificare se la luminescenza emessa dal cristallo sotto l’azione 
della radiazione y, fosse effettivamente responsabile dell’effetto osservato, ho 
ritenuto opportuno, come prima prova, fotografare lo spettro della luce emessa. 
tid è stato fatto usando un Hilger E2. Lo spettrogramma riprodotto nel testo ri- 
vela una banda compresa fra 4000 e 5000 A, con un massimo d’intensitä intorno 
a 4460 À. Essa & la più intensa delle tre bande rivelate fotograficamente da 
L. Ror (?) usando un cristallo di antracene e radiazioni di alta energia. 


(°) L. Rorx: Phys. Rev., 75, 983 (1949). Cfr. pure F. B. HARRISON e G. REYNOLD: 
Phys. Rev., 79, 732 (1950). 
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Per ottenere il mio spettrogramma, che è del resto assai poco intenso, sono 
state necessarie ben 72 ore di esposizione, mentre con il FM si ottengono imme- 
diatamente effetti notevoli. 

Ulteriori prove analoghe a quelle di AGENO e collaboratori (1) eseguite per 
accertare che l’azione registrata dal FM & dovuta eclusivamente alla lu- 
minescenza sono: 

1) togliendo il cristallo o la vaschetta con la soluzione e lasciando davanti 
alla finestra del FM la sola sorgente dei raggi y non si è notato alcun aumento 
della corrente di zero. Cid dimostra che i raggi y non hanno azione diretta sul 
fotocatodo; 

2) interponendo fra il cristallo o la vaschetta irradiati, e la nenn del FM 
un cartoncino nero la corrente anodica si è ridotta alla corrente di zero. Cid 
dimostra che il contributo degli elettroni estratti dal cristallo o dalla solu- 
zione per azione dei raggi y incidenti è addirittura trascurabile; 

3) sostituendo al cristallo o alla vaschetta con la soluzione una lastrina 
di piombo di spessore tale da emettere sotto l’azione dei raggi y un numero 
di elettroni Compton maggiore di quello emesso dal cristallo o dalla soluzione 
in esame, non si è notato aleun aumento della corrente di zero. 

Assunta arbitrariamente come unità di misura per la corrente di Jumine- 
scenza quella fornita dal cristallo di antracene alla tensione di 900 V fra gli 
estremi A e B (fig. 1), sono state eseguite misure comparative su una massa 
di g 7,520 di soluzione in xilolo rispettivamente di naftalina e acenaftene alla 
concentrazione di 0,18 g/em?. Successivamente si sono effettuate misure sulle 
quantita di xilolo puro e delle sostanze organiche che entrano come compo- 
nenti in 7,520 g di soluzione. 

Le misure relative ai liquidi e alle sostanze organiche eseguite in condi- 
zioni geometriche identiche sono riportate nella seguente tabella. 


TABELLA I. 


| | 
| 
| 


Corrente 


SOSTANZE | Grammi ; ; 
diluminescenza 


3,760 | l 


EIS divantracene aaa ee] 3 

Soluzione di naftalina in Lin : 7,520 0,107 

Xilolo puro . ; 5,972 0,012 | 

Naftalina in sta 1,548 0,066 | 

Soluzione di Acenaftene in xilolo . 7,520 0,138 | 

Xilolo puro 6,000 0,013 
1,520 0,071 


. Acenaftene in el 
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Da essa risulta che tanto lo xilolo quanto la naftalina e l’acenaftene, irra- 
diati separatamente, presentano una debolissima luminescenza, la quale cresce 
quando essi sono uniti in so- 

17108 luzione. 

Si è anche esplorato il com- 
portamento della intensità della 
luminescenza delle soluzioni al 
variare della concentrazione. 

Il grafico della fig. 5 rappre- 
senta l’intensitä della corrente di 
luminescenza in funzione della 
concentrazione. 

Come si vede, una piccolissima quantita di acenaftene è sufficiente per 
produrre un netto aumento della luminescenza (praticamente nulla per lo xilolo 
puro); questa aumenta rapi- 
damente al crescere della con- 
centrazione, raggiunge un 
massimo e poi lentamente de- 


0 0,04 008 012 016 0,18 
Concentrazione grammi/em? 


Fig. 5. - Soluzione di acenaftene in xilolo. 


cresce per concentrazioni via 
via crescenti (10). 

E stato anche esaminato 
il comportamento di una so- 
luzione di naftalina in xilolo 
(0,18 g/em*) al fine di poter 
stabilire un confronto diretto 
tra le misure di altri autori 
e quelle qui riportate in fig. 6. 
L’andamento della nostra cur- 
va à lievemente diverso da 
quello trovato da altri au- 
tori (11), Appare infatti, dalla 
fig. 6, che per concentrazioni 
maggiori di 1,5 % l’intensità 
della Juminescenza rimane 
praticamente costante, mentre 


quella degli autori menzionati Fig. 6- a) soluzione di naftalina in xilolo; b) par- 
ha un andamento crescente. ticolare della curva. 


0% 257 10 MIS 20) 25 35 40 50 


Concentrazione mg/en® 


(1°) Alla concentrazione 0,18 g/cm?, la soluzione di acenaftene o di naftalina & prati- 
camente satura. 

(4) M. AGENo e collaboratori: loc. cit (1); H. KALLMANN e M. Furst: Phys. Rev., 
81, 853 (1951). 
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Dei fatti osservati si possono dare, a mio avviso, due spiegazioni, se pure 
qualitative. 

La prima pud ottenersi, penso, seguendo l’analogia con certi fenomeni di 
assorbimento, nella regione del primo ultrarosso, da parte di varie sostanze 
disciolte in tetracloruro di carbonio (1?) (1%). Si pud allora pensare che la lumi- 
nescenza provenga quasi esclusivamente dalle molecole disciolte e si presenti 
cosi intensamente nelle soluzioni, perch i movimenti delle molecole, di qua- 
lunque natura esse siano, si effettuano pit liberamente, essendo qneste mag- 
siormente allontanate. 

La seconda spiegazione possibile, analoga a qnella suggerita da H. KALL- 
MANN e Furst (5), si pud basare su un’analogia con i fenomeni di fluorescenza 
sensibilizzata (Woop, RASETTI), pensando che la luminescenza sia prevalen- 
temente prodotta dalle molecole del solvente, ma si manifesti intensamente 
solo quando, per l’aggiunta del soluto, le molecole del solvente consegnano a 
quelle in esso disciolte, la loro energia di eccitazione. 

Nop sembra che al momento attuale si disponga di elementi sufficienti per 
accreditare l’una o l’altra spiegazione 0, eventualmente, qualche altra diversa 
da esse. 

Sono allo studio soluzioni di altre sostanze organiche in xilolo, benzolo 
ed altri liquidi. 


Ringrazio il prof. CONVERSI per le utili discussioni avute sull’argomento. 


(2) L. BERTI: Nuovo Cimento, 5, 44 (1948). 
(13) A. Crocoxe: Nuovo Cimento, 5, 489 (1948); 6, 500 (1949); 7, 373 (1950). 


SUMMARY (*) 


The paper reports experiences performed in order to study the luminescence of 
solutions due to +-Trays, and the intensity variation of luminescence with changing 
concentration. Two possible qualitative explanations of the luminescence are reported. 


(*) Kditor’s translation. 
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Summary. — The Argon isotopie composition has been measured in many 
italian fumeroles and soffioni. A remarkable enrichment of A! has been 
found, certainly due to the K! decay in the inner of the earth. With 
the help of these data, the problem of the origin of atmospheric Argon 
is considered. 


Introduction. 


Natural gases have been studied for a long while to get informations on 
the inner of the earth. Chemical analyses of natural gases (volcanic gases, 
fumeroles, drill-hole gases and so on) have been made by many authors (1), 
but as far as we know the isotopic composition has not yet been considered. 

Rare gases are always present in some amount in natural gases. Attention 
must be stressed to rare gases as invaluable sources of information on consti- 
tution, origin and evolution of the earth atmosphere. 

Argon (?) isotopes are very interesting in this respect. The K“ disinte- 
gration was postulated by VON WEIZSÄCKER (?) to explain the great amount 


(1) For a review of these data see: U. SBORGI: Ann. di Chim. Appl., 82, 395 (1942). 
(2) The Argon content in the atmosphere is 0,93% in volume. Its isotopie com- 
position, as given by A. O. NIEr: Phys. Rev., 77, 789 (1950), is 
AS 0.337 + 0.001 % 
ASE 0.063 + 0.001 % 
AX 99.600 +- 0.001 % 
(3) C. F. von WEIZSÄCKER: Phys. Zeits., 38, 623 (1937). 
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of A‘ in the atmosphere. As a matter of fact K* (which is contained in 
natural Potassium as 0.0119 %) has been shown later to decay not only by f 
disintegration, but also (*) by K-capture, the final product ofthis second process 
being just A. Furthermore Potassium content in the earth has been shown (°) 
to be in such a large amount to justify the A“ content in the atmosphere. 

A question still open is the mechanism by which A“ accumulated in 
the atmosphere (%). There is still some doubt whether a great amount of the 
actual A“ was already contained in the primitive atmosphere, or whether it 
was continuously evolving in the geological time. In this respect the A4/A35 
ratio in natural gases is a source of information for the A% amount which 
continually goes in the atmosphere. 


Experimental part. 


We starked our work in volcanie gases. It is a difficult task to capture 
true volcanic gases avoiding air contamination. A few attempts have been 
made to catch gases from Etna, Stromboli and Vulcano, introducing a vacuum 
tight tank as deep as possible in the vulcano gaps. A device was used to 
open and to close the tank right inside the gap, in order to prevent any air 
contamination. In all collected samples the enrichment factor 


N= [A20/A%6] Perea Ao! | Sophos : 


was so close to one to believe the collected smokes were not really coming 
from the inside of the volcano, but probably were some eireulating air in 
the gaps. 

For this reason we left voleanic gases for fumeroles and drill holes where 
the capture is much easier and safer (expecially for the experimentalist). As 
two samples provided by the interest of the Larderello Corporation were 
found (7) to have a 7 value larger than one, we collected further samples from | 
the most important italian fumeroles and « soffioni » ($). 


(4) F. C. Tmompson and $. Rowranps: Nature, 152, 103 (1943). 

(5) H. E. TAreL: Journ. Geoph. Res., 55, 329 (1950); F. BIRCH: Journ. Geoph. Res. 
56, 107 (1951). 

(6), J. H. J. Pooun, C. F. G. Decaney: Nature, 167, 680 (1951); A. ROSTAGNT: 
Geof. pura e applicata, 18, 128 (1950). 

(7) G. Boato, G. CARERI, G. Nencını and M. SANTANGELO: Ann. di Geofisica, 4, 
111 (1951). 

(8) « Soffioni » are high temperature and pressure steam wells, somewhat like geysers. 
Some of them carries to the atmosphere a hundred tons of water vapour per hour. 
For details see: P. Grnort Conti: Journ. Chem. Edue., 4, 281 (1927); U. SBORGI: Atti 
della Reale Acc. d'Italia, 5, 667 (1934). 
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«Soffioni» samples have been collected and purified by the Larderello 
Corporation. Fumerole gases have been collected by us with à bell immersed 
in the bubbling water. Water and carbon dioxide, the most abundant com- 
ponents (°), were both adsorbed on the spot. Rare gases were extracted in 
the laboratory in the usual way. A quantity of nearly 100 em? at 20 mmHg, 
suitable for isotopes analysis, was obtained working on one cubic meter of 
original gases at standard conditions. 

A Nier-type mass spectrometer built in this laboratory (1°) was used for 
isotopes analysis. For every sample the full mass spectrum in the range 35 


x 240 230 2.20 210 200 


Fig. 1. — The typical mass spectrum of atmospheric Argon, and the background (dotted 
curve) in the mass range 36 to 40. It is a plot of the ion current in arbitrary units 
versus the magnet current in Ampere. 


to 41 has been measured and compared with an atmospheric Argon sample 
used as a standard. In this mass range the background had a constant value 
after the copper analyser tube had been heated at 200°C for a few days. 
The background has been always measured before and after a sample has 
been introduced, and subtracted in the proper way. Due to the high in- 


(?) Except for water, always present as vapour or also in condensed form, carbon 
dioxyde constitutes 90-99% of these gases. No appreciable oxygen quantities are 
generally found, instead Nitrogen and rare gases are always present. Other gases, 
like Hydrogen, Hydrogen sulfide, Methane, are sometimes present in small amount 
as reported in reference (+). 

(20) G. CARERI and G. NENCINI: Nuovo Cimento, 7, 64 (1950). 
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tensity of mass 38 background peak, the Argon measurements have been 
limited to the 36 and 40 peaks. Magnetic scanning was used with an 
1500 V accelerated ion beam. Molecular flow was surely established, an 
important condition to make quantitative analysis without a calibration of 
the spectrometer. As a proof of the good performance of our instrument, we 
obtained for the A*/A%* ratio in atmospheric Argon 294 + 3, to be compared 
with the Nier’s 296 +1, obtained with a carefully calibrated instrument. 
A typical mass spectrum of ours is reported in fig. 1. 

Results of our analyses are reported in table I. The quoted error is mainly 
due to the reproducibility, as for this purpose some samples have been measured 
again after a few months. 


TABLE I. 
= "i 7 = | 
SOURCE | [A20/A39] gas/ [AP / AP] atın. | 
FUMEROLES 
Agnano (Napoli) 1.15 (a) 
Pozzuoli (NapolDereme ae | 1.12 
Bullicame (Viterbo) ram. vr 1.00 
Montemiccioli (Volterra) . . . . . . | 1.99 | 
Rapolano (Siena) EEE 0 1.26 | 
| | 
| SOFFIONI 
| Larderello N. 1 (6) (Pisa)... . . | 1.22 
| Larderello » 2 Da, cl 1.19 
| Larderello Seo 5 | o>, rt 1.16 
| Larderello » 4 DT ee 1.24 
| Larderello » 5 5. CANON 1.28 
| Castelnuovo N. 1 DIR SA RE 1.12 
Castelnuovo » 2 D pelle re 1.40 
Serrazzano N. 1 et ER 1.60 | 
Serrazzano » 2 De ue [era 
Travale MP es ie 1.21 
|  Sasso ‘Souk Pale 1.25 
Lagoni Rossi CR te 1.45 
Lago Tassinaie DR Re 1.25 
| Monterotondo ee ER Be 1.45 


(a) Experimental error on the enrichment factor + 4%. 
(b) A progressive number has been used to identify different neighbor sources. 
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Discussion of the results. 


With a glance at table I, one sees the observed enrichment to be sometimes 
very remarkable, compared with the usual variations of isotopic composition 
of elements in the earth. 

The A*/A% ratio in all samples is never less than the atmospheric one, 
and this seems quite reasonable from the following considerations. The atmo- 
sphere of to-day is generally assumed to be completely of secondary origin (11), 
namely the result of the earth outgassing from the time when gravitational 
field and termal conditions were such to prevent gas excaping far from our 
planet. Although some outgassing still remains to-day, and volcanic pheno- 
mena are an example of it, the most important ‘quantity of gases has been 
developed probably during the earth’s solidification. 

Let us consider Argon isotopes A and A®. Due to the same chemical 
properties and the light mass difference, we can think they have had the 
same story in the process of earth outgassing. Also without any guess about 
the A#/A® value in the atmosphere at the beginning of the earth, one must 
say, a8 a result of K* decay into A%, that ratio has always been less in the 
atmosphere than in the inside of the earth. In other words, the A*/A* ratio 
in the atmosphere has alway been rising, but with some delay with respect to 
the inside of the earth. 

Obviously one must always consider only the average value of the ratio 
A*/A8° in the inside of the earth. In spite of our results, it can also happen 
that some natural gases have a 7 value less than one, and this is due to the non 
homogeneous Potassium distribution. The figures of table I show remarkable 
variations of 7 depending on the gas source. 

We believe this is due to the non uniform geological constitution of the 
surface crust, to the different path of the gas coming to the surface and also 
to some air contamination inside the earth. It is difficult to day to explain 
these variation without a knowledge of the gas path, and also of data con- 
cerning » of the same source in different times. 

To answer the question about the mechanism of A“ accumulation in the 
atmosphere, one should extend this kind of analysis to all the earth’s surface, 
to know more about the quantity of natural gases excaping up to day. It is 
however interesting to point out that the Soffioni of Larderello supply about 
4.4 : 105 grams of A* per year, to be compared with 6.6 - 101? grams (12), which 


(1!) H. JEFFREYS: The Earth (Cambridge, 1929), p. 312; H. Brown: The Atmo- 
spheres of the Harth and Planets (Chicago, 1950), p. 260. 
(2) F. A. Panrtu: Neue Wege exact. Naturerk., pag. 1 (1939). 
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is the total A% content in the atmosphere. Although it is not possible to 
think that such an intense outgassing happens all aroud the earth, the above 
figures allow us to believe that a remarkable amount of A% is still flowing in 
the atmosphere. 
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RIASSUNTO 


La composizione isotopica dell Argon & stata misurata in un gran numero di fuma- 
role e di soffioni italiani. Si é trovato. un notevole arricchimento in A, sicuramente 
dovuto al decadimento del K* contenuto nell’interno della Terra. Alla luce di queste 
misure viene considerato il problema della genesi dell’ Argon atmosferico. 
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Summary. — Centres of symmetry, mirror planes, and rotation axes can be 
recognized in suitable cases through a statistical study of the intensities 
of the X-ray reflexions. The number of essentially distinct space groups 
that can be determined by X-ray methods is thus raised from about 60 
to 215, out of a total of 219. Statistical methods may also be of use in 
the study of pseudosymmetry. Applications to ephedrine and haemo- 
globin are described briefly. 


Certain symmetry elements in crystals (lattice centring, glide planes, screw 
axes) may be readily recognized by the systematic absences they produce 
among the X-ray reflexions. Others (centrosymmetry, mirror planes, rotation 
axes) have no such obvious effect on the X-ray reflexions, and until recently 
were thought to be undetectable by X-ray methods. For crystals with poorly 
developed externa] forms, therefore, there was often considerable doubt about 
the proper assignment of crystal class and space group. Recently, however, 
it has been shown that a statistical study of the intensities of the X-ray 
reflexions will often reveal the presence of the second type of symmetry element, 
and the number of crystals of uncertain symmetry can be correspondingly re- 
duced. 

If the contributions made to the structure factor by the various atoms. 
of a non-centrosymmetrie structure are plotted on an Argand diagram an 


(*) A paper given at the Second International Congress of Crystallography, Stock- 
holm, 1951. 
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irregular polygon is formed, bearing a striking resemblance to illustrations 


of Brownian motion (fig. 1). The line représenting the resultant structure 
factor can have any length from zero up to a maximum equal to the sum of 


90° 


Fig. 1. — Contributions of individual atoms added vectorially to form the structure 

factor of a non-centrosymmetric crystal. The broken line is closely analogous to 

illustrations of Brownian motion. (The device shown is a draughting machine modified 
for the rapid and convenient graphical caleulation of structure factors (1).) 


the atomic scattering factors. Except for crystals with few atoms in the unit cell 
or for very low values of sin 6, however, the probability of the structure factor 
coming anywhere near this maximum is very small. In various ways (?) (?), 
one being the analogy with Brownian motion just mentioned, it can be shown 
that the average value of the intensity of reflexion, J = F?, is equal to the 
sum of the squares of the atomic scattering factors, &, and that the proba- 
bility of I lying between J and-1 -+d/ is 


(1) P(Dal = S-exp(— 1/2] dl. 


Integration (4) shows that the fraction of the intensities less than or equal 
to I is given by . 


(2) N(I) = 1 — exp 2/2]. 
() D.-C. Pnrczrps: Journ. Sei. Instrwm., in course of publication. 
(2) A. J. C. Wırson: Nature, 150, 152 (1942). 
. (3) A. J. C. Wırson: Research, 2, 246 (1949); Acta Oryst., 2, 318 (1949). 
(4) E. R. Howeııs, D. C. Parezrps and D. ROGERS: Research, 2, 338 (1949); Acta 


Oryst., 3, 210 (1950). 
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The variance, the mean square deviation of J from its mean value 2, is easily 
seen to be %?, and in fact a more exact expression can be derived (5), free 
from restriction on the number of atoms in the cell: 


V= 22227, 
where 2, is the sum of the fourth powers of the atomic scattering factors. 
For centrosymmetric crystals the situation is rather different, as the atoms 
occur in pairs whose combined contribution to the structure factor is restricted 
to the real axis. The mean value of J still remains 2, but by analogy with 
one-dimensional Brownian motion or otherwise it can be shown (?) that the 
probability distribution of the intensities is now given by 


(4) IUD) CUE == 7 D) ER pie PS OWE. 
the fraction of intensities less than J is 


(5) IND) = Grae (L/P) 


2 


where erfx is the error function, and the variance is 22°, or more exactly 
(6) V=227— 3d). 


Expressions (4), (5) and (6) are quite different from (1), (2) and (3) in their 
mathematical form, and experience with over a dozen substances has shown 
that distinct indications of the presence. or absence of a centre of symmetry 
can be obtained by their intelligent use. The qualification intelligent is im- 
portant, because some judgement is required in deciding (i) whether the sub- 
stance is sufficiently complex, and (ii) whether a sufficient number of reflexions 
has been estimated to make the application of statistical methods justifiable. 

One substance to which the methods have been applied is the alkaloid 
ephedrine. The fact that the mean intensity of reflexion is equal to & may 
first be used to place the relative estimated intensities on an absolute scale. 
Since thermal vibration reduces the intensities of reflexion by the factor 
exp [— 2B sin? 6/A?], a graph of log (<I)/2) against sin? 0 will be a straight 
line, whose intercept on the ordinate axis gives the factor required to convert 
the estimated relative intensities to the absolute scale, and whose slope gives 
the effective value of the temperature parameter B. (The effective value 
may be increased above that for temperature alone by crystal imperfections.) 


(5) A. J. C. Wırson: Research, 4, 141 (1951). 
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Three such graphs are shown in fig. 2; they are linear within the experimental 
error and determine the absolute scale within à few percent. Fig. 3 shows 
N(z) as a function of z (where z = 7/2) for the same two zones, together with 
the theoretical curves from equation 3 (lower curve) and 6 (upper curve) 
There is no doubt that the 010 projection of this substance is centrosymmetric 


[001) (-Eph HCI 


[010] 
(-Eph HBr 


01 0-2 03 04 05 06 07 08 09 


Fig. 2. — Graph of log (Z/(1)) as a function of sin? 9 for two zones of ephedrine: hydro- 

chloride and one zone of ephedrine hydrobromide. The plots are closely linear, and 

from their intercepts the factors required to place the arbitrary visual estimates of 

intensity on an absolute scale can be found. Their slopes give the effective Debye 
temperature factors. 


+. 
60 we) rat 


N(2)(%) 


Fig. 3. — Fraction (in %) of reflexions with intensity less than 2% of the average 

value. The curves are theoretical, the upper one for a centrosymmetric crystal or 

projection, the lower one for a non-centrosymmetrie crystal or projection. The 

experimental points are for two zones of ephedrine hydrochloride, and show that the 

010 zone corresponds to a centrosymmetrie projection, and the 001 zone to a non-centro- 
symmetrie projection. The space group P2, is thus confirmed. 
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and that the 001 projection is non-centrosymmetric, confirming the space 
group P2,. / : 

The differences in intensity distribution may be applied either to the 
general reflexion hkl for the detection of real centres of symmetry, or, a8 has 
just been described, to zones. of reflexions for the determination of centro- 
symmetric projections. In projections, however, the statement that the ave- 
rage intensity is equal to 2 has certain exceptions (°), and these make possible 
the statistical determination of symmetry elements that cause atoms to coin- 
cide in projection. For example, in projection on to a mirror plane the atoms 
coincide in pairs, and for the corresponding zone of reflexions only the crystal 
must be regarded as containing half the actual number of atoms, each having 
a doubled atomic scattering factor. Since the scattering factors are squared 
in 2, the average intensity in this zone is doubled. In general, when the atoms 
in projection coincide in groups of n, the average intensity in the corresponding 
zone or row is multiplied by n. Comparison of zone or row averages with each 
other or with the general average thus makes possible the detection of a 
considerable number of symmetry elements. The use of statistical methods 
and systematic absences together make possible the unambiguous determination 
of 215 of the 219 effectively distinct space groups (7), whereas the latter alone 
determine unambiguously only about 60. 

An interesting possibility not yet fully explored is the study of pseudo- 
symmetry. If a structure is pseudocentrosymmetrie in the atomic positions 
one would expect the distribution curve to change from the type of equation (5) 
to the type of equation (2), and the variance to decrease from 22? toward 2, 
as batches of reflexions of successively larger mean sind are examined. 
(Other types of pseudocentrosymmetry might produce the opposite variation). 
Through the kindness of Mr. E. R. HOWELLS and Dr. M. F. Prrurz I am 
able to reproduce two N(z) distribution curves for the 100 projection of horse 
methaemoglobin (fig. 4). The upper curve shows the theoretical distribution 
for a centrosymmetric crystal, and the experimental points for the low-order 
reflexions are in good agreement with it. The lower curve is the theoretical 
distribution for a non-centrosymmetric crystal, and the experimental points 
for the higher-order reflexions are in fair agreement with it, though there is 
some irregularity. These curves confirm the conjecture, based on other evi- 
dence, that this projection of methaemoglobin is pseudocentrosymmetric, 
though it cannot actually be centrosymmetric. 

The limited time allowed for this paper makes it impossible to discuss the 
practical methods of carrying out the statistical investigations. I shall there- 


6) A. J. C. Witson: Research, 3, 48 (1950); Acta Oryst., 3, 258 (1950). 
(") D. Rogers: Research, 2, 342 (1949); Acta Cryst., 3, 455 (1950). 
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fore close by summarizing: in many cases the statistical methods make pos- 
sible the conversion of relative intensity data to the absolute scale, make 


0 20 40 60 80 100 


Fig. 4. — N(z) curves for horse methaemoglobin, 100 projection. The crosses are from 

low-angle reflexions (spacings greater than 6 A), and the circles for higher-angle 

(spacings less than 6 A). The curves are theoretical. The projection is thus actually 

non-centrosymmetric, in accordance with the space group, but would appear centro- 
symmetric if examined at low resolution. 


possible the detection of symmetry elements not revealed by systematic 
absences, and may give indications of pseudo-symmetry (not necessarily pseudo- 
centrosymmetry). Knowledge of the absolute scale and pseudosymmetry can 
be very helpful in the early stages of structure determination. 


RIASSUNTO (*) 


I centri di simmetria, i piani di simmetria, e gli assi di rotazione Si possono Spesso 
identificare con lo studio statistico delle intensità delle riflessioni dei raggi X. Il numero 
dei gruppi spaziali nettamente distinti che puö essere determinato con l’impiego dei 
raggi X si eleva cosi da eirca 60 a 215 su un totale di 219. I metodi statistiei possono 
anche essere impiegati nello studio della pseudosimmetria. Si deserivono brevemente 
aleune applicazioni del metodo all’efedrina e all’emoglobina. 


(*) Traduzione a cura della Redazione. 
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Cosmic Ray Evidence on the u-meson Magnetic Moment. 


D, C. PEASLEE (*) 


Washington University, St. Lowis, Missouri (**) 


(ricevuto il 1° Dicembre 1951) 


Summary. — A survey is made to determine what information is provided 
by the intensity-depth relation of cosmie rays underground on the pos- 
sibly anomalous character of the |1-meson magnetic moment. Energy 
losses due to ionization, bremsstrahlung and electron pair creation by 
normal and anomalous Dirac particles are discussed. It is concluded 
that if the magnetic moment in terms of the Bohr magneton is (1 + à), 
the only limitation provided by present experience at depths of 3000 mw.e. 
or less is |ö| S 1. 


1. — Introduction. 


The intensity versus depth relation of cosmic rays underground can be 
interpreted solely in terms of electromagnetic interaction of the u-meson (t), 
down to depths of 3000 meters water equivalent. Agreement with obser- 
vation is obtained by assuming energy losses appropriate to a normal Dirac 
particle. The question arises whether this agreement is sufficient evidence that 
the u-meson has a normal magnetic moment, since particles of anomalous 
moment have exceptional radiative properties at high energy. A negative 
answer is indicated by the following argument. According to present views, 
only particles of compound structure have anomalous moments. This com- 
pound structure may involve one or more (virtual) particles with a maximum 
associated Compton wavelength 2,. The anomalous moment effects will then 
be cut off for processes involving momentum transfers Ap = h/A, in the rest 


(*) Now at Columbia University, New York. 
(**) Assisted by the joint program of the ONR and AEC. 
() S. Hayakawa and S. TOMONAGA: Prog. Theor. Phys., 4, 287 (1949). 
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system of the original particle. For example, on this account the anomalous 
moment terms are of little interest with regard to divergences in quantum 
electrodynamics, since they will have a negligible effect unless it is assumed 
that there exists a truly elementary particle with anomalous moment; but such 
a particle would represent almost as grave a philosophical difficulty as the 
divergences themselves. 


2. — Ionization losses. 


The cross section for energy loss by electron knock-on of a normal u-meson 
goes as In #,; the anomalous moment introduces a term proportional to E,, 
which should dominate at high energies except for the cutoff mentioned above. 
If the magnetic moment is (1 + 6), with u, the normal moment, the cross 
section (?) for knock-on with transfer of energy eE, to the electron in the 
laboratory system is, for Æ#, > uc’, | 


(1) o(E,,e)de = dar, (l/s —1)de, 


where r, = e?/ue? is the «classical radius» of the meson. The integrated 
cross section for energy loss is clearly proportional to E, if no cutoff is im- 
posed. For the cutoff the longest Compton wavelength involved is presumably 


that of the u-meson; thus in the center of gravity system 
(2) 12 Apluc = (2p/ue) sin (1/2)6 , 


where the momentum p in the ¢.g. system is simply related to E, and m, the 
electron rest mass. The cutoff (2) may be inserted in (1) by means of the re- 
lation & sin? (1/2)0; the integrated cross section for energy loss then becomes 


D (E öl 2) 2 | B, x, = 0.8 Lee |m , 
. yj = 0° (a )ocr, y v 
(3) o(E,) zur, lutetim  E,— 00. 


These asymptotie forms indicate a maximum energy loss per unit path lengt 
of | DE/Dd« | <1.9-10-? MeV em2/g. Of course the normal ionization loss is 
<10, which is the only 


present in any case and will be dominant for | 0 
limitation imposed by experimental agreement with intensity-depth curves 
:aleulated for 6 — 0. 

The ionization produced by à fast neutron can also be crudely estimated 
in this way. Here the cutoff wavelength is presumably that of a r-meson;: 


€ 


€) H. C. CORBEN and J. ScuwinGer: Phys. Rev., 58, 953 (1940). 
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the maximum energy loss occurs at Hy ~ 0.5[m,M/mlce? — 1.3 - 10° GeV and 


is about 5:10 MeV cm?/g, less than 1% of conventional «minimum joniz- 
Ation ». 


8. — Bremsstrahlung. 


The anomalous bremsstrahlung can be treated by the well known method 
of virtual quanta (5). The leading term in the anomalous cross section is 


By "max 
(4) o(E,,, e)de = dta(Zr,) (1 — E) de | dÆ, | dr/r , 
euc?/1—e) r 


min 


with ~ —1/137. The integrals in (4) describe the Compton scattering in the 
meson rest system of a virtual photon of energy Æ, from the nuclear Coulomb 
field; the lower energy limit is obtained from transformations relating the 
initial and final Compton scattering energies, and the meson rest system with 
the laboratory system. The radial integral is characteristic of the Weizsacker- 
Williams method; when no special restrictions are imposed max 4,H,/E,, 
where 2, is the Compton wavelength of the meson. The lower limit rin is 
roughly equal to the cutoff (nuclear) radius R of the Coulomb field or 2,, 
whichever is larger; for the u-meson R = 2,Z1% to a good approximation, 
while in the electron case A>R for all Z. Thus the integral over r is 
In(gE,/E,2%3), where g is a number of order unity; without cutoff, the integral 
over E, is of order Æ,, so that the integrated cross section for energy Joss is 
of order Æ*. | 

If the eutoff wavelength for the anomalous terms is taken to be that of 
the u-meson, the energy integral is to be extended only over E, = ue? and 
becomes 


az ( /1—(3/)2 « 
Eu Zils 10 < 2 13 
(5) | dE, In (gE,/E,213) | 1 ji (gEul ue?) E<24)9 ; 
euc?/(1—e) 0 eZ 2/3 s 


This cutoff therefore causes a reduction to normal type bremsstrahlung; and 
further reduetion to the form In (const) for very high E, is brougth about 
in the usual way by electron shielding of the nuclear Coulomb field, which 
euts it off at r,,,— @271%, with a, the Bohr radius. The integrated cross 


max 
(2) ©. F. v. WEIZSÄCKER: Zeits. f. Phys., 88, 612 (1934). 
(4) S. B. Barporr and R. THomas: Phys. Rev., 59, 621 (1941). 
(5) J. L. POwELL: Phys. Rev., 75, 32 (1949). 
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section for energy loss by anomalous bremsstrahlung is thus of order 
o(E,) ~ 0.120to(7r,)2E, In (gF) 


x E,|Z"3uc? 
{ 6 F == | u J 
N (3.8: 10/2285; 


where the smaller value of F is always to be chosen. The corresponding 
cross section for normal bremsstrahlung, always present in addition, is 
a(Zr,)°E, In gF. Because of the 0* dependence, the anomalous bremsstrahlung 
is a somewhat more sensitive index than the ionization loss; in view of the 
agreement with 6 =-0 down to 3000 mw.e., | ö| can certainly not be as 
large as 10, but a value of order 1 is not excluded. Intensity measurements 
at depths >10° mw.e. would in principle determine the bremsstrahlung loss 
more precisely and might reduce the upper limit on | 6 |. 


4. — Pair production. 


A similar discussion applies to the process of electron pair creation, which 
for mesons is comparable in importance to bremsstrahlung and would consi- 
derably exceed bremsstrahlung for a particle as heavy as the proton. An 
estimate of the normal cross section for energy loss by pair creation (the ano- 
malous cross section can again be neglected for |6| <1) is given by the 
Weizsäcker-Williams method for BE, > uc?: 


"max 
(7) Opair (2/7) a | 6, dE, | dr/r, 


"min 


where a, is the cross section for electron pair creation by a photon of energy E,. 
For high E,, the range of E, greatly exceeds 183 Z-3 me?, so that the constant 
high-energy value can be used for o, throughout. In this case rn is deter- 
mined by the longest Compton wavelength involved, which is that of the 
electron; and "ax is A„B,/E,. Integrating E, only over the range for which 


Tmax = "min; one has 
(8) Te (8/z)(xZr,)’(m/u)gE, In(183Z-¥8) , 


where r, is the electron classical radius and m its mass. Taking g—1, the 
ratio of pair production to bremsstrahlung at high energies is 


= ~~ (œ/x)(u/m) . 


mro\ In (183771) 
| In (3.8 - 104 Z72/®) 


(9) = Opair/Obrems GS (me 


ce 8 
UT u) 
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The dependence of R on (w/m) shows that pair creation is negligible in 
comparison with bremsstrahlung for a primary electron, is of comparable 
importance for a meson, and is probably the dominant effect for a proton. 
With R~ 0.5 for a u-meson, the energy loss per unit path length due to both 
pair production and bremsstrablung is about 2:10 E, em?/g, which is to 
be compared with the roughly constant ionization loss of about 2 MeV cem?/g. 
Therefore only the effects of normal ionization loss need be considered for 
mesons with | 6| S1 at underground depths of 10? mw.e. or less. 


RIASSUNTO(®*) 


Si compie un’indagine per determinare quali informazioni si possano trarre dalla 
relazione intensita-profondita dei raggi cosmici sotto terra sul carattere eventual- 
mente anomalo del momento magnetico dei mesoni y. Si diseutono le perdite 
d’energia dovute alla ionizzazione, alla bremsstrahlung, e alla creazione di coppie di 
elettroni da parte di particelle di Dirac normali e anomale. Si conclude che, se il mo- 
mento magnetico espresso in magnetoni di Bohr & w,(1 + 6), la sola limitazione for- 
nita dalle presenti esperienze a profondita di 3000 my,o 0 minori, è |d| S 1. 


(*) Traduzione a cura della Redazione. 
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Summary. — Cosmic rays underground are treated as entirely caused by 
the electromagnetic interactions of y-mesons. Simple approximations 
are given for the underground spectra of u-mesons, virtual photons, 
soft showers, and nucleons. Comparison with experimental data is satis- 
factory on the absolute rate of slow neutron and nuclear star production, 
and on the charge ratio of star primaries, if the charged primaries are 
at least half protons. Disagreement is found with the observed variation 
of star production with depth. Comparison with observation suggests 
that the cross section for nuclear excitation by y-rays decreases on 
the average more slowly than 1/H from its large measured value 
o,/A = 1-10-28 em? per nucleon at E = 200+ 300 MeV. This mechanism 
therefore predominates over meson-proton knock-on in production of 
the underground nucleonic component. Measurement of slow neutron or 
star production over large differences in depth would yield more detailed 
information on the behavior of o,/4 at high energies. Measurement of 
the size-frequeney distribution of soft showers at high energies would 
yield information about the production of neutral 7-mesons. | 


4. — Introduction. 


Of the cosmic radiation at sea level, only u-mesons and neutrinos can pene- 
trate to depths of about 20 meters water equivalent or more. There are 
accordingly only three sources of underground cosmic rays: i) electromagnetic 
interactions of y-mesons; ii) non-electromagnetic interactions of u-mesons; 
iii) neutrinos. Only the first of these is of any importance, and can be used 


(*) Now at Columbia University, New York. 
(**) Assisted by the joint program of the ONR and AEC. 
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exclusively to interpret the intensity-depth curve down to 3.000 mw.e. (4). 
If this is assumed as the minimum mean free path for a non-electromagnetic 
interaction, the corresponding cross section is & S 5:10-% em? per nucleon... 
Direct measurements (2) have indicated o = 2.3 10% cm? for fast mesons;. 
and the value calculated from the specifically nuclear potential (3) necessary 
to explain the capture of slow negative mesons is a ~ 107% (p/uc)? em?, where p 
is the meson momentum and w its rest mass. The same formula in the same 
units holds for excitation by neutrino absorption, considered as an inverse 
p-decay process, where p is now the neutrino momentum, 

Thus only the electromagnetic interactions of y-mesons will be considered 
below; and the accompanying note indicates that it will be sufficient to con- 
sider the meson as a normal Dirac particle, even if this is not quite the case. 


2. — Underground meson spectrum. 


The vertical cosmic ray intensity as a function of depth shows a consi- 
derable change in slope ad about 300 mw.e.: at shallower depths it varies (4) 
roughly as d-*5, at greater depths as d7?5. This has been interpreted (5%} 
as reflecting a competition between nuclear absorption and u-decay for the 
parent 7-meson high in the atmosphere. In the 7-meson rest system, the 
u-meson kinetic energy is negligible in comparison with its rest mass; there- 
fore iu the laboratory system Æ,Jue— E,fm.e. Accordingly the relative 
probability of r-u. decay is 
(1) P(E) ={ (m_¢?/H,,)1/te} lk of + (m,c?/E.\1]rce; = Bi(E,+.B), 


u 


B = huc/ot , 


where À is the mean free path for absorption by air nuclei of a r-meson of 
energy E_ and mean life + in its rest frame, and o is the air density at the ave- 
rage height of production of u-mesons. 

If the x-mesons are originally distributed according to a power law E,’"!, 
the integrated u-meson intensity will vary with depth as d”’ far above and 
d-’~* far below the critical depth corresponding to the energy loss B. This 
statement assumes a linear energy loss with depth by the u-mesons, which 


(1) S. Hayakawa and S. Tomonaca: Prog. Theor. Phys., 4, 287 (1949). 
(2) E. Amazpr and G. FIDECARO: Nuovo Cimento, 7, 535 (1950). 
(8) J. A. WHEELER: Rev. Mod. Phys., 21, 133 (1949). 
(4) V. C. Wırson: Phys. Rev., 58, 337 (1938). 

(6) K. I. GREISEN: Phys. Rev., 73, 521 (1948). 

(8) S. Hayakawa: Prog. Theor. Phys., 3, 199 (1949). : 
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is justified because the critical depth is considerably less than 10? mw.e., 
where bremsstrahlung losses first begin to become important. The average 
value obtained from intensity-depth measurements is y = 1.65, which is in 
accord with the observation (7) that the high energy differential spectrum of 
vertical mesons at sea level can be fitted by (H, + Æ,) 2%, where H, ~ 20 uc? 
is roughly the ionization loss by mesons arriving from the top of the atmo- 


2 - 
sphere. It is also in agreement with the value y —1 3 recently deduced (°) 


for the primary radiation at high energies. 
For depths at which linear energy loss is valid, the differential vertical 
spectrum underground is 
[| n.(E,)dE, = mP(E, + od)dE,[E, + od)" 
(2) oo yelb 
Ny = 2.0(ue?)*° cm? 87! sterad-!, 


where o ~ 2ue?/m w.e. and d is the depth measured from the top of the atmo- 
sphere. The low energy mesons at sea level that deviate from this formula. 
will all have been absorbed at less than 10 mw.e. underground, and decay 
in flight is negligible. The coefficient x, is evaluated by fitting (2) to the sea 
level measurements (°) at Æ, ==10uc?, where P~1. To discuss the integral 
spectra, crude approximations will be used throughout to preserve a 
perspicuous analytic form at the expense of precision, which is not to be 
expected to a high degree in any case. The total vertical intensity obtained 
by integrating (2) is 
ee 
(3) 1,4) = | n,(B,) AH, = nl (d/y(od) , 


0 


where it is clear that J, > 1 as od/B — 0, and J, > y/[(y 4- 1)od] as od/B — =. 
The simplest form with these properties is 


(3a) J,(d) ~ PI + 1/y)od], 


which will be assumed as a lowest order approximation (in excess). 
Under the form (3a) the depth D of the «knee » in the intensity — deptb 
curve is related to the constant B by 


(4) Bw Dl + 1/y)> 108 we, 


Janossy: Cosmic Rays (Oxford, 1948), p. 171. 
V. NEHER: Phys. Rev., 83, 649 (1951). 
Rossi: Rev. Mod. Phys., 20, 537 (1948). 


(7) L. 
() H. 
(°) B. 
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using D~ 300 mw.e. Detailed calculations (1%1!) show that at angle 6 to 
the vertical 

(5) B(0) ~ B/eos 6. 


The source of local secondary radiation underground is the total meson 
spectrum integrated over all angles 9 of deviation trom the vertical. Assuming 
only ionization (linear) energy loss at depths d < 10% mw.e., one has 


(6) I,(d) = I,(d/cos 0) = [no/y(od)’][(cos 0)’J,(d/cos 6)] , 


where by (3a) and (5) J, is actually independent of 0, and the integral of (6) 
over all angles becomes 


(7) I(d) = 2an,P[(1 + 1/y)od]/y(y + 1)(od) . 


It is also desirable to have an expression for the differential spectrum n(E,) 
integrated over all angles; a simple approximation of suitable general form is 


270) dE. 
[Eu + (y + Doar 
which has the properties that it integrates to (7) and agrees asymptotically with 


the direct angular integration of (2) in the limiting cases Æ, > 0, co. With 
the insertion of numerical values given above, these become 


(8) n(E,) dE, = P[( + 1/y)od] , 


(Ta) I(d) = 0.92 d-258[1 + 2.7 10-3 d] 1 em= 8-1, 


and the normalized form 


if 
(8a) (de ee 


7 Bu) dL, = 13.7 d*[B, + 3.6 dy * dB, , 


where Z,, is in units of wc? and din mw.e. Equations (7a) and (8a) are the 
basis of the remaining discussion. 


8. — Secondary soft radiation. 


Soft showers are initiated by electron knock-on and bremsstrahlung; 
u-meson decay is negligible in the dense medium underground. As indicated 
in the accompanying note, electron pair creation by w-mesons is similar in 
character to bremsstrahlung and about half as large in magnitude. The cor- 


(2) K. GREISEN: Phys. Rev., 76, 1718 (1949). 
(41) K. GREISEN: private communication. 
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responding cross sections per atom for energy loss are in simplest approximation 


Or(B) AE = 2ar)Z(me?) dE/ EB? 
Op(E)dE =- 4a(Zr,)Æ, In (3.8 + 1047-2) dE/E 
op(E)dE = 0.50,(#) dE 


o(E, E,) = 6g + Op + Op, 


where E, E, are the shower and meson energies, respectively, and 7, and r, 
are the classical radii of electron and meson. In gp there is cancellation of 
two factors of 2 from the fact that two showers are initiated, cach on the ave- 
rage with half the energy of an equivalent bremsstrahlung shower. For 
Showers of energy E the production rate R,(E) in g1s-1 is given by 


fos) 
Es {Caco,] _ 
(10) [RABITAIAE >, te 
p Br 71—.027E) 
= (N/A) | o(E, E,) dE una Hess | -10=dB/E®, 
E 


with # and d in we? and mw.e. The two sets of coefficients in the brackets 
correspond to CaCO, and Pb, respectively, and will be carried throughout as 
typical of light and heavy materials. If the measured size-frequeney distri- 
bution should show any pronounced excess above these curves, it would 
presumably constitute evidence for a strong nuclear source of soft showers, 
such as the decay of neutral 7-mesons. This excess would be apparent only 
for showers of energy E > m,c?. 

Since the cross section for nuclear excitation by y-rays appears to be 
rather high, it is desirable to have an expression for the average photon spectrum 
accompanying an underground meson. The photons are both virtual and 
real, arising respectively from the Coulomb field of the meson and the soft 
showers. The virtual photon spectrum accompanying the average meson is 


CO 


9 de 
(41) y(E,)dE, = = «(dE,/E,) [ In (E,/E,)i(E,) dE, = = «(dE,/E,L)(E,). 


2 
The integral L(E,) occurs repeatedly, and the following forms are found nume- 
rieally: 


0.86 In (3.6d/E,) E,< 3.6d 


(12) L,(E) tal + E,j3.6d)-"® E,>3.6d, 


a - Il Nuovo Cimento. 
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where the low- and high-energy expressions do not coincide in the intermediate 
region H,~3.6d. The logarithmic term is in error by less than 10% for 
.002 < H /d< 2 and suffices for most purposes; the second formula is needed 
only for rare events of very high energy. For the real photons the approxi- 
mation is used that a shower initiated by an event of energy FE’ has a path 
length ('?) for y-rays of energy E, that is X(E,)dE, = 0.57X,E'dE,/E? where 
X, is the radiation length of the material traversed, X," = 4«(Zr,)2N,/A x 
xIn (18321). The corresponding spectrum of real photons per average 
meson is 


co Eu 
(13) y’(B,) dB, 0.57(X,N,/A)(dB,/H%) | i(E,)dE, | B'o(E', B,) aH! = 
By Ey 
TO | 


~ (dE, JE?) | L(E,) + 1.4: 1042.59 — y-155)| 


l1.9-10-2J 
y = (1 + #/3.6d) , 


where as usual the brackets refer to CaCo, and Pb, and Æ,, d are measured 
in uc? and mw.e. In evaluating (13), the rough approximation is made that 
the average value of E, in the integral of type L is En 3.6d. 

Because of the different dependence on #,, the virtual photons are do- 
minant for energies E, = 5 ue? 0.5 GeV, the real photons at lower energies. 
For nuclear excitation, proton knock-on makes some contribution, so that 
the average distribution of secondary electrons per meson might also be com- 
puted. This distribution would be given by 0.76 y"(H;) and for energies 
Ls =, ue? would be only a few percent so intense as the primary meson com- 
ponent; since the energy region above 100 MeV is the most interesting for 
consideration of nuclear excitation, the secondary electron component will be 
neglected throughout. 


4. — Nuclear excitation. 


Nuclear excitation of energy E by proton knock-on occurs at a rate 


R,(E) g-4s-1, where 
: co 


(4) -R,(B) AB /I(d) = 2m N,Z/A)(m/u)(m|M)(4B/E?) | i(2,) dB, = 


E 


3. 
hs 10-"(dE/E2)[1 + 43.641". 


13-2) 


(2) B. Rossı and K: GREISEN: Rev. Mod. Phys., 13, 240 (1941). 
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This rate represents à minimum estimate, since it can ouly be increased by 
the knock-on of virtual m-mesons (#); since the 7-mesons are present only 
a fraction of the time, however, equation (14) should probably be correct in 
order of magnitude. The excitation rate by y-rays is RE) gs", with 


(15) R,(E)dE/I(d) = Noy(E)o,(E)/4 , 


where y(E) =y’(E) + y"(B), and o,(E)/A is the effective cross section per 
nucleon for y-excitation. 

Recent measurements (14) indicate that in the range #, ~ 200+300 MeV 
the total cross section for excitation is of order N,o,/A — 10-4 g-1 em?, a re- 
markably high value. At these energies comparison of (14) and (15) shows 
that A, exceeds R, by an order of magnitude. Hence measurement of size- 
frequency distribution of nucleon cascades underground would yield inform- 
ation about the behavior of o, at energies above the limits of presents machines; 
such measurements would in principle extend up to energies where o,(E) starts 
to fall off faster than 1/#. 

As an index of nuclear excitation one may consider the average number 
of slow neutrons m, produced per unit length (in g/em?) when a meson tra- 
verses an underground absorber. This number is 


co 


(16) md) =| [Y(R, + R)I@NEIE)AE, 


E 


$s 


where H,~ 0.1 uc? is the average of a «slow » neutron relative to its bound 
state, and the yield Y is the fraction of initial excitation energy that is even- 
tually converted into slow neutrons. Certainly Y<1, and Y = 0.1 seems 
probable as a lower limit. Two possible types of variation will be assumed 


for o,: 

(17a) 10-* 
N,0,(2)/A = 

(175) 10-4(2/E) , 


in both cases H>2, and o,(E) =0 for H<2. Measurement of m, as a 
function of d would be evidence as to which of equations (17) is closer to the 
truth. 

For comparison with observation it is convenient to define o,(d)/A, an 
effective cross section per nucleon of slow neutron production; and for pro- 


(3) C. FRANZINETTI: Nuovo Cimento, 7, 384 (1950). 
(4) S. Kixucur: Phys. Rev., 81, 1060 (1951); Bull. Am. Phys. Soc., 26 (3), 7 (1951). 
(5) R. D. Mrzcer: Phys. Rev., 82, 260 (1951). : 


68 D.:C. PEASLEE 
duction in Pb at d = 30 mw.e., 


(184) 82.Y © 10529 sem? 
g,(30)/A = m,(30)/No = 
(18b) Wea ea IDES cern’ 


A preliminary measured value (1) can be titted by (18a) with Y — 0.5, but 
(18b) would require Y~2, which is impossible. This is an indication that 
at high energies the excess of k, over R, becomes even greater than the order 
of magnitude found for excitation around 200-300 MeV. The underground 
nucleonie component is therefore primarily due to y-excitation of nuclei and 
only secondarily to preton knock-on. 


5. — Stars in emulsion underground. 


Nuclear stars in emulsion have been measured underground (17) at depths 
d = 33, 44, 70 mw.e. from the top of the atmosphere. Very little trend with 
depth is evident, and the formula 


(19) Ry. (d) = 2.0 : 10-7 4-57 g-1 g-1 


represents about the maximum decrease with depth consistent with the stated 
errors. The average of all the measurements is 


(20) R, (55) ~1.5- 10-8 gs. 


About 4% of the primaries appear to be charged. The observed stoppings of 
slow 7-mesons cannot be readily interpreted, as it is not clear to what extent 
the slow rs are locally produced in the material immediately surrounding 
the emulsion. A r-meson with a residual range of 100 u of emulsion has a 
90% chance of decay in a path length of 3 meters. It therefore seems likely 
that the observed rate is strongly influenced by apparatus geometry. 

The energy required for a star of 3 or more prongs is taken (18) as 
H=1Aypc?, The stars are produced directly by y-mesons (knock-on and 
virtual quanta), by soft showers (real quanta), and by nucleons in locally 
produced cascades. For the nucleon cascades, a crude diftusion treatment 
is developed in the Appendix, which results in the following approximate 


(6) R. D. SARD: private communication. 
(7) E. P. GEORGE and J. Evans: Proc. Phys. Soc., 63, 1248 (1950). 
(8) R. Brown et al.: Phil. Mag., 40, 862 (1949). 
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path lengths for neutrons and protons of energy E initiated by an event of 
energy Es: 


2 A„(E, Ey) dE ~ 0.42 4,E,(dE/E®){1 + 3E,/(E, + 1.6E)] 
1,(E, Ey) dE ~ 0.42 7,E,AE/[EIE 4- E,)], 


where the ionization loss by a proton per geometric mean free path A, is 
1.4 pe? [ CaCO, | 
E el for 1. The average number of nucleons of ener 
a0 | Pb | ne 5 er 
accompanying a meson is then 


ve) 


(22) N,,9(H) AB — dB | { An oH, Ey)[Rx(Eo) + R,(e)]/L(d) } ¢ 
E 


In the emulsion the rate of star production is somewhat less than the rate 
of nuclear excitation, since a nucleus with sufficient energy to produce a star 
of 3 or more prongs does not always do so. For an order-of-magnitude esti- 
mate it is assumed that the average rate of visible stars is 1/2 the rate of 
excitation to energies = 1.4 wc, regardless of the mode of excitation. Then 
the star rates with meson, y-ray, neutron and proton primaries are 


R,, = (1/2) | [R.(B) + R,(B)IdE 
1.4 
(23) R,, = (1/2) | R,(E)dE 
1.4 
== 
Run, (1/2)1(d)(No0,/-A) | N,,(E)dE, 
1.4 


where N,0,/A — 1.0 -10-? em?/g is the average geometrical cross section of 
emulsion. 

In evaluating R,,, in (23), it is sufficient to regard the factors external 
to 2,, as relatively slowly varying with E, so that the correction terms 
involving Æ, can be approximated by averaging over /,,,(2, E,) alone, for 

<H< oo. With E, =1.4 ve, this yields 


AE, Es) dE ~ 1.68 : 0.42 7, E,AE/E® 


(24) i 
1(E, E,) dE ~ 0.69 : 0.42 A,B d B/E? 


so that 30% of the star producing nucleonic component consists of protons. 
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This extimate, although only order of magnitude, indicates that the proton 
component is not to be ignored. With this approximation, expressions (23) 
are evaluated under assumptions (17a) and (175), resulting in table I. The 
fractions of the total star rate due to each Source are given, as well as the 
traction of charged primaries f, = f, + f 


1 ee 


TABLE I. 
Variation 
R,,(55) 3070 N ik Im 1» fe 
mw.e. 
Approx. (17a) ae Er d--87 .09 .09 .58 .24 .33 
Approx. (17b) |" 0.7.1078 ers ta 19 SG | 19 38 


Comparison with the observations (19), (20) leads to the following con- 
clusions. The absolute star production rate suggests that the truth is inter- 
mediate to approximations (17a) and (17b); that is, the relatively large value 
observed for o,/A at H = 200 —300 MeV is maintained over a considerable 
energy range, falling off more slowly than 1/#. This has the corollary that 
meson-proton knock-on is of generally negligible importance as an interme- 
diate mechanism in production of underground cosmic rays. The fraction 
of charged primaries is relatively insensitive to assumptions and is about 35%, 
of which at least half appear to be protons. The variation in depth cannot 
be made to ‘agree with the very small experimental value, and an assumption 
intermediate to (17a) and (175) would yield something like d-1 as the most 
probable calculated variation. 

The variation in depth would be more significant if it could be determined 
at greater depths, such as shown in table II. The values of f, are also included 
and are seen to show very little variation. 


TABLE II. 
R,(a) f(a) R,(b) f.(b) 

| = _— m = = ——— = meet uns ne na! Seel pr 
| d — 100 m w.e. la a a yee .31 SONORE Er pr) 
| | 
| 5 PRE a Seale 
| 150 1.0 + 10-8 o-1 5-1 31 1.7 - 109 g-1 8-1 33 

200 Se NORE ote 30 ea I oot i Opec eset ses 31 


309 | AG T0 ag SA Fr Ole ait ‚64: 10207871 29 2 
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6. — Summary and conclusions. 


Cosmic rays underground are considered as due to the electromagnetic 
interaction of u-mesons. Tractable approximate expressions are obtained for 
i) the flux of underground u-mesons; ii) size-frequency distribution of soft 
showers; iii) rate of nuclear excitation and production of slow neutrons; 
iv) production of nucleon cascades and stars in emulsion. 

The dominant feature in nuclear effects is the large cross section o,/A for 
+-excitation of the nucleus, which has been measured to be of order 1 - 10-25 em? 
per nucleon at HE = 200 to 300 MeV. Comparison with measurements on 
slow neutrons and star production underground indicates that o,/A falls off 
more slowly than 1/E for higher energies and hence remains relatively large 
over a wide range. Other modes of nuclear excitation, such as meson-proton 
knock-on become negligible in comparison. 

The agreement calculation with observation is satisfactory except for the 
depth dependence of nuclear stars; this difficulty would be illuminated by 
measurements at greater depths than hitherto attempted, for which calculated 
values are given. Measurement of nuclear excitation by means of stars or 
slow neutrons over wide depth intervals would yield more definite information 
about the variation of o,/A at high energies. 

Measurement of the size-frequency distribution, of soft showers at energies 
of a few hundred MeV and above would also provide information on the cross 
section for neutral z-meson production at these energies. In all such appli- 
cations, the underground y-mesons are being utilized essentially as a naturally 
provided betatron of low intensity but unlimited energy range. 


The author takes pleasure in thanking Professor R. D. SARD for an intro- 
duction to these problems and for many interesting and helpful discussions. 


APPENDIX 


In the absence of ionization loss, the diffusion equation for a nucleon 
cascade reads 


co 


(Aly © dN(E, x)jdx =— N(E, x) + | GE, E')N(E',x)dE', 
E 


where N is the nucleon density and x the distance in units of y,, the geometric 
mean free path. For the reproduction funetion a convenient phenomeno- 
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logical form is 
(A2) G(E, E)dE'—= (n + 2)(n + 1)(E'— EyE’"1dE', 


where the normalization insures energy conservation in a nuclear collision 
process with E’ the incident energy and E the energy parameter of the 
emitted particles. Here n is the single adjustable parameter, since there is 
only one experimental fact to be reproduced; namely, that the mean free 
path for absorption of the atmospheric N-component is 24,. The Mellin trans- 
form of (A1) is 


(n + 2)! s! 
(n—+s-1)! 


pas N(s, x) = [—1 + 2C(s)IN(s, x) , 


| DON (s, 2)/dr = 


(A3) 


N(s, 2) = [(E/B,)'N(B, 2) ak . 


0 


For an initial power spectrum N(E,0) = E-* with s — 1.65, the value 
n — 2.5 makes the bracketed coefficient in (A3) become — 1/2 in accordance 
with observation. 

For energies at which ionization loss by the protons becomes )important, 
the equations become 


DN ,(E, @)/d0 = — N,(E, #)+-[ GB, B'AG/2N,(E", 2) + (1/2)N,(B", #)) a" 


(A4) oo 
3N,(B, 2)[d%0 = —N, (EB, ») + [G(B, BAG/2IN(E", 0)+(1/2)N,(E', #)]dE'+ 
+ E.DN,„(E,z)/dE, 


where Æ, is the (assumed constant) ionization loss per mean free path. The 
factor of 1/2 inside the integrals assumes equal probability for neutron and 
proton emission in nuclear collisions. Since these equations are difficult to 
solve, consider the easier set obtained by inserting an additional factor (E/E,) 
in the ionization loss term. This will cause a considerable underestimate of 
ionization losses but is at least more realistic than (A3). The transformed 
equations are then 


| dW, (8, z)/ de = — N,(s, 2) + C(s)[N,(s, x) + N,(s, x)], 
(A5) dN,(s, æ)/dæ EE N,(s, 2) SF C(s)[N,(8, x) ci N,(8, x) | zz 
| — (8 + 1)(E.JEo)N,(8, x) . 


The solution with N,(s, 0) + N,(s, 0) =1 corresponding to an initiating event 
of energy E,, which reduces to the solution of (A3) if H, — 0, is 
| N„(8, æ) | se { C/(u + 1) | 


Hates) idee aa 


(A6) 
| u(s) = O(s) —1— 6/2 + V0? + 62/4, d=(s+1)EJB,=(s +1)6)- 
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The total path length at energy Æ is the inverse Mellin transform of 


co 


| N(s, x) dx and can be approximated by the residue at 89, which has the alge- 
0 

braically largest real part of any pole position. For 6 = 0, this pole of 
l/u(s) is at so —1; and s, = 0.95 even for 6, — 10, so that s, ~1 can be 
taken as a good approximation throughout. The residue at s, is considerably 
more sensitive to 0); but simple numerical approximations good to DOS 
the range 6, — 0 to 10 are 


(AT) | AE, Ho) AE = 0.42 4,[1 + 304/(1 + 1.66,)] 2, AE/E? , 


A,(E, By) dE — [0.42 ,1/(1 + d)]E, dE/E®, 


where the geometrical path length A, is introduced explicitly. 

The path lengths (A7) are based on an approximation that reduces the 
actual ionization loss by a factor H(B)/E,<1, where E(E) is the average 
energy of a nucleon cascade in passing from energy E, to E. The crudest 
estimate is # — E; in higher approximation this will be corrected by factors 
of order In (E,/E), ete., which are secondary in order of magnitude. Using 
just the lowest approximation, one has 6, > 6,H,/E, with the simple result 


a | An(E, Ho) dE 0.42 1,[1 + 3E./(E. + 1.6E)]E,dB/R: , 


AE, Ey) dE — 0.42 4,E,dE/E(E + E,). 


RIASSUNTO (*) 


I raggi cosmici nel sottosuolo sono trattati come dovuti eselusivamente alle. intera- 
zioni elettromagnetiche dei mesoni u. Si danno delle semplici approssimazioni per gli 
spettri sotterranei dei mesoni u, dei fotoni virtuali, degli sciami molli e dei nueleoni. 
Il confronto coi dati sperimentali è soddisfacente per quanto riguarda la misura asso- 
luta della produzione dei mesoni lenti e delle stelle nucleari, e il rapporto di carica dei 
primari producenti stelle, se i primari carichi elettricamente sono almeno per meta 
protoni. Non c’è accordo con la dipendenza osservata della produzione di stelle dalla 
profondita. Il confronto con l’osservazione fa pensare che la sezione d’urto per l’ecci- 
tazione nucleare da parte dei raggi y decresca in media più lentamente di 1/Æ dal suo 
notevole valore misurato o,/A — 1: 10-28 em? per nucleone a Æ = 200300 MeV. Questo 
meccanismo predomina quindi, nella produzione della componente nucleonica, nel sotto- 
suolo su quello di knock-on mesone-protone. Misure di neutroni lenti o di produzione 
di stelle a grandi differenze di profondita fornirebbero maggiori chiarimenti sul com- 
portamento di o,/A alle alte energie. Misure della distribuzione ampiezza-frequenza 
degli sciami molli alle alte energie fornirebbero dati sulla produzione di mesoni 7 neutri. 


(*) Traduzione a cura della Redazione. 
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Summary. — This paper develops the theory of an improved type of resistor 
network, for the numerical calculation of an harmonie function (in two 
dimensions, or in three with axial symmetry) with given values on the 
boundary. The advantage of this type of apparatus on those hitherto 
described is that it takes into account terms of an higher order in h 
(h = side of the mesh), so permitting a much better accuracy with the 
same value of h. It is shown that, also for the usual type of network, 
the reasoning based on the analogy with the electrolytic tank is mis- 
leading, and that the correct method of calculation leads, in the case 
of axial symmetry, to values of the resistances. which, in the vicinity of 
the axis, differ by an appreciable amount from those adopted by preceding 
Authors. A method is indicated for imposing the boundary condition 
on a curved line not passing through the junctions of the network. 


1. — Introduction. 


As it is well known, the problem of solving numerically Laplace’s equation 
AV — 0 with given boundary conditions is of paramount importance in many 
fields of physics and engineering. Since an analytical solution can be found 
only in a minority of cases, three methods have been developed for an appro- 
ximate numerical determination of the function V: the relaxation method (*), 
the electrolytic tank and the resistor network. The first is a method of nume- 


(7) R. V. SOUTHWELL: Relaxation Methods in Engineering Science (Oxford, 1940); 
Relaxation Methods in Theoretical Physics (Oxford, 1946). See also: G. LIEBMANN: 
Adv. in Hlectronics, 2, 101 (1950), where a neat account can be found also of the 
methods of the electrolytic tank and of the resistor network. 
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rical computation whose approximation can be pushed as far as desired, but 
at the cost of a considerable sum of work and time. The last two methods, 
based on electrical analogies, once constructed the necessary apparatus enable 
to solve a variety of problems with little labor and accuracy sufficient in most 
practical cases. It has been pointed out recently (2) that the resistor network 
can give results much more accurate than an eleetrolytie tank and presents 
many advantages over the latter. 

The resistor network so far considered is a net of electrical resistors di- 
sposed according to the pattern of 
fig. 1: they have all the same value 
if the equation in cartesian coor- 
dinates 


ee eo, 


must be solved: they have instead 
values graded according to à cer- 
tain law if the equation to be sol- 
ved is 


that is, the equation satisfied by an 

harmonie function of the cylindrical 

coordinates 7, 2, 9, in the case of Fig. 1. 
axial symmetry (V independent 

of 9). This is, for instance, the equation of the electric (or magnetic) potential in 
an electron lens. It can be shown that the electric potential » at any junction 
of the net (which can be measured by means of a potentiometer) equals the 
value of the function V at the corresponding point of the field, provided that 
the fourth order derivatives of V, in the case of (1), or the third order ones, 
in the case of (2), can be neglected. 

In this paper we will show that the accuracy can be much improved by 
adding to the network appropriate resistors along the diagonals, as shown in 
fig. 2. In this type of network the lowest derivatives that are disregarded are 
those of eighth order in the case of (1) (cartesian coordinates) and of the fourth 
order (but with a very small coefficient) in the case of (2) (cylindrical coordi- 
nates). It is found that, in the case of (1), the diagonal resistances must be 
four times larger than the others; in the case of (2), the law is more compli- 
cated, and, besides, also the non-diagonal resistances must be recaleulated by 
means of recurrence formulae. But these calculations, made once and for all, 
enable to construct any resistance network of this type: these results are 
collected in a numerical table (tab. I). 

Of course the increase in accuracy obtained in the network with diagonals 
would be deceptive if the boundary conditions (which we take of the Dirichlet 


AS | 


(2) G. Lippmann: Brit. Journ. Appl. Phys., 1, 92 (1950). 
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type V= V,) were not realised with a corresponding exacteness. Since the 
boundary line, along which must be V=V,, in general will not pass exactly 
through the junctions, a method is indicated, more accurate than that hitherto 


employed, by which the junctions 
nearest to the boundary will assume 
automatically potential values diffe- 
ring from V, by the appropriate 
small amounts (n. 6). 

We have based our calculations 
on Taylor developments, avoiding 
any resort to the analogy with the 
electrolytic tank which, in our opi- 
nion, is misleading. In fact, if our 
method is applied to the ordinary 
network (without diagonals) for cy- 
lindrical coordinates, it is found 
(n. 4) that the resistances must be 
graded according to a law differing 
from that adopted by de PACKH (*) 
and by LIEBMANN (?): the difference 
is negligible at points far from the 
axis 7 — 0, but not so in the vici- 
nity of the latter. For the resistan- 
ces along the axis, we find a de- 


finite value intsead of the infinite value given by de Packh’s formula. 
Our values of the resistances for the network without diagonals for eylin- 


drical coordinates are also collected 
im table I. 


2. — Diagonal network for cartesian 
coordinates. 


In this section we suppose tha- 
the unknown function V(x, y) sa- 
tisfies equation (1) (cartesian coort 
dinates). We immagine drawn in 
the wy plane a lattice of square 
meshes of side h, whose corners cor- 
respond to the junctions of the re- 
sistor network. Let P, (fig. 3) re- 
present any of these junctions (not 
contiguous to the boundary). For 
the sake of simplicity, we have sup- 
pressed from this drawing as weli as 
from the following ones all the re- 


sistors except those connected to the junction P,, and have represented these 


(8) D. C. pe PACKH: Rev. Sci. Instr., 18, 798 (1947). 
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1 
—] 


with straigth lines instead of the customary symbol. The resistances along 
the horizontal and the vertical sides will have all the same value 1/C, those 
along the diagonals the value 1/H (so C ad H mean the conductances: they 
are the same throughout the network). The potentials of the junctions 
Ps, Pı, Pi,..-, will be denoted by v,, v,, v;,.... and the values of the V fun- 
ction at the corresponding points of the xy plane will be called V,, V,, Wy.... 
The differences 0, — V,, v, — V,, etc., will be the «errors ». It is impossible 
to give a general expression for the errors, but an insight on the accuracy of 
the method can be gained by evaluating the error & —%,— V, at P,, when 
the errors at all eight neighbouring junctions are assumed to be zero. 
Kirchhofl’s law applied to the junction P, gives 


(3) ey nn 0) Hm 05 Tu Fu en), 
and therefore, if v, — V,, vj = Vi, etc, % — Vo + & 
(4) A CV; + Koss %»+V,— 4V,) EHE, EN) 


4(C + H) 


Now V,, Vi, ete., can be expressed as a Taylor series by means of V and its 
derivatives at P,: denoting by h,, h, the coordinates of a point with respect 
to P,, the series can be written 


co d d n 
5 RN Wr 
ö) ne Sn, +R), 


(where it is understood that the derivatives are to be taken at P,). In de- 
veloping this formula we must take into account that, because of (1), 


2) > \? 
> ni <a 


(if applied to V or any of its derivatives) and therefore, for any integer k, 


0 ar coy” 


Applying (5) to the eight points for which h,, h, have the values 0, + h, 
and substituting in (4), we find that many terms cancel out and we are left 
with 


co ham dam V 
VA Ar VG ae V, SF Ve 4Vo = 4 = (4m)! ain ? 


ham dt" Va 


y; ae Ve Ae Va 1 = 4V, == ti In ( ayn (4m)! yam 0 
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Therefore, if we disregard all derivatives higher than the eighth, formula (4) 
reduces to 


af ht } d4 V hs 8 V 


8 — - = ) Se ER er 
(8) École 


From this expression it is seen that we can get rid OF the fourth order part 
of the error by making 


(9) Cate 


that is, giving the diagonal resistances a value four times that of the others. 
The principal part of the error then becomes 

hs dV Ge OP 
~- 10080 da® 10080 2% ' 


(10) 


This value must be compared with the error in a network of the usual type, 
without diagonals. We can find the latter by putting H—0 in (8), which 
gives for the principal part of the error, without diagonals, 


ht DEV ht DV 
24 dat 24. dyt 


(11) — 


It may appear that, with the cost of the materials used for the diagonal 
resistors, one could instead enlarge the usual type network and so adopt a 
smaller value of h, but it is readily seen that h would be reduced only by a 
factor V2, and therefore e by a factor 4. The advantage gained by reducing 
the error from the fourth to the eighth order is in general much more consi- 
derable (except, of course, in regions where V happens to be very well re- 
presented by a polynomial of the third degree). Besides, it is to be expected 
that the doubling of the number of resistors connected to each junction con- 
tributes to a better smoothing out of the random errors caused by the resistors. 
tolerances. 


8. — Diagonal network for cylindrical coordinates. 


We turn now to equation (2) and proceed along the same line of the pre- 
ceding section, except that the values of the resistances of a mesh depend 
on the distance of it from the (horizontal) axis r—0. So the eight resi- 
stors connected to any junction P, (with r 40) are of five different values, 
as shown in fig. 4, where the letters C, G, G', H, H' denote, as before; 
the conductances. The values of each of these quantities is a function of 7, 
that is of the order number n=r/h of the row. The row n —0 will be 
considered later (n. 5): for the moment we suppose n 4 0. 

Kirchhoff’s equation (3) for the junction P, is now replaced by 


OO ri) G(V, — Vo) 4 G'(v; Vo) + 


+ Ho; + u — 20) + Ho, + vf —20)) = 0 
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and if ww — V, + €, v, =N, ete:, the error is given by 
I ” 1 1 ce 0 


(12) | RS BT EEE 
2C+6+@+2H+M) 


where we have called / the quantity 


Ber OP on) Lan M) GV. = V,) 
Vs chet 2) (Vian 2 Vo) 
To calculate this expression we develop V,— Vy, Vi— V,, ete., with for- 


~ 


mula (5) (where « and y must now be replaced by z and r) limited to the fourth 


order terms. Then the numerator I of (12) comes out as a rather cumbersome 
expression containing the eight derivatives (at P,) 


(13) Ves Ver , Ver , Ver I Ve ’ Vea , Vs; ’ Vie , 


dr 
lated by equation (2) and those which can be deduced from it by differen- 


tiation, namely 


[where = (27) ; ete.) These quantities are not independent, being re- 
Po 


VW, = Vie 
Ve Re Ve == Ar Ar V,[r? 
(14) RE ve vr 2 V,,fr? — 2V,|r® 


Var = ar = 2 om Vr/r? = V.[r® ‘ 
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By means of these formulae we can eliminate from the expression of J all 
derivatives with respect to z, and we are left with the following formula, where 
the derivatives are all independent (and n = r/h) 


| il il jt il it ul 
I 2 eis > = "2 
1-6 G + ( as ina) E | or. 5 siekes + seal H I 
LÉO AS 
=p = F a) c Br 
DT IE 1 h? 
ñ | ’{ = t Op > 2 
(15) a ß 1 G T | N =e = dal SF É I Fi) H | al c\ 2 ! Vr+ 
2 2 Ch: 
Fo e—_@—(4+ a + [4 \#4 et 
n N, n| 3! 


ht 
GG Ge 20] eee 


Now, if we suppose already determined the conductances @', H’ of the (n—1)-th 
row, we can determine C, G, H (conductances of the n-th row) so as to 
eliminate the first, second and third part of J (and therefore of e). In fact, 
putting equal to zero the coefficients of V,, V,,, Vir, we get the following 
three linear equations in the three unknowns C, G, H: 


(16) 12134 2 (24n3 — 12m? = 12m AA (4 272— 1) ¢ = 
— 12n3G'+ (24n3 + 12n? + 12n + 11)H’, 


(17) 6n2G — (2n— 11) (nr ED 6 C2 nr UE 
(18) nG— (4n + 2)H + C =n@'— (4n—2)H’. 


We start the solution by eliminating C. It is convenient at this point to 
introduce the following abbreviations: 


L, = 24n? + 36n?—16n+ 9, P, = 247 —n, 
M, = 48n? + 24n? + 16n—9, N, =12n?+ 6n?+ n. 
Then the equations for G and H take the form 
PG eo DB = PG! 11 ’ 
N,G— M,H =— N_,G'+ M_,H', 
and their solution is 
=a (LnN-, Ss PM.) = Dur ar L,M_,)H 
L,M,—P,M, : 
Tr. PAN ar Wn Goat (PMN he)! 
LINE ag P,M, ; 


(19) GS 


(20) = 
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Once calculated G and H, formula (18) gives 


(21) O =n@— n@ — (4n—2)H'+ (4n + 2H. 


These formulae enable to calculate, step by step, all the conductances of 
the network when those of the first row are given. Inn. 5 we shall show how 
the first row of « vertical » and «diagonal » conductances (that is, G and H 
for n = 0) can be calculated once the value of the conductances along the 
axis (C for n = 0) has been arbitrarily chosen. In this way, using formulae 
(35), (36), (19), (20), (21), we have computed the values (to five significant 
figures) of all the resistances of the first 26 rows, assuming for those along 
the axis the arbitrary value 1/0, — 100000 ohm: the results are collected 
in table I. Of course all resistances could be multiplied by an arbitrary 
constant. 


TABLE I. — Values of the resistances. 
With diagonals | Without diagonals 
| n ] T j — 
| 1/C | 1/G | 1/H 1/C | 1/G 
| | | | 
tu ro 100 000 | 14706 71,438 | 100000.0 | 25000.0 
El 6 578.9 4 612.5 18 248 | 15625.0 10135 
(hyo 3 321.5 27134, So 10 714 | 7865.3 6214.5 
PRE 2221500100 RO lala! 7 620.8 5245.9 | 4466.9 
Veet 16673 | 1491.4 5 921.1 3 934.8 3483.4 
1295 1 334.2 | 1 217.9 4 843.3 3147.9 | 2853.9 
RT LL ON aa 1029.30 17740982...) 26232 | 24167 
7 953.01 891.40 | © 3552.0 2 248.5 2 095.5 
8 ace: || 786-1 UNE rs 134-49 0 3) 4 11:967.4 | 113496 
9 TATA E 703.09 2 804.8 1748.8 | 1655.3 
10 667.02 | 635.94 | 2537.8 1573.9 1497.9 
11 606.31 | 580.51 | 2317.3 1 430.9 1 367.8 
12 555.80 533.97 2 132.1 1 311.6 1 258.5 
PAS 513.03 494.35 1 974.3 1210.7 1 165.3 
ee 476.34 460.20 1 838.2 1 124.3 1 085.0 
MIE 444.60 430.47 1 719.7 1 049.3 1 015.1 
| 16 416.78 404.35 1 615.5 983.74 953.64 
heey 392.24 381.21 1523.3 925.88 899.19 
| 18 370.49 360.59 1441.0 874.43 850.62 
19 351.01 342.07 1 367.1 828.50 807.03 
| 20 333.47 325.37 1 300.5 787.02 767.64 
| 317.61 310.22 | 1 240.0 749.51 731.95 
22 303.18 296.42 1 184.9 715.46 699.44 
23 290.01 283.80 1 134.5 684.36 669.70 
24 277.93 272.20 1 088.2 655.86 642.34 
25 266.83 261.52 1 045.6 629.60 617.16 
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If the resistances are chosen according to this law, in formula (15) all terms 
disappear except those of the fourth order and the error becomes 


ht o*y 


20+ 44 G'—8(H+H!) 
FOOTAGE GRIT EH): 


(23) & 


The values of « for a few values of n are given in table II. It will be noted 
that they are very small and tend to zero when n increases, as it is to be 
expected because, very far from the axis, conditions approach those of the 
cartesian coordinates -(n. 2). 


TABLE II. 
| 
I 1 | 2 5 10 15 | 20 
| _ — — — 
| | 
| @ = 0.0529 | 0.0016 | —0.0034 | —0.0015 | — 0.0009 | — 0.0006 


4. — The network for cylindrical coordinates without diagonals. 


Now we shall avail ourselves of the calculations of the preceding section 
for deducing the theory of the network without diagonals in a better approx- 
imation than that adopted by de PACKH and LIEBMANN. 

Formula (15), in the case of a network without diagonals, that is with 
H= H'= 0, becomes 


he ee Ls il 3 a : | il RTE | 
es) 1=|6—0 + ( = 4 sn c RV,+ a+ |a+ @ (2 - ss) Cle: ie 
GC ES } a 
== G— G+ Ale Vere + [20 + Oe Ge Were 
n|3! 4! 


If @' is given, we can determine C and @ so as to suppress the first and second 
order terms. This leads to the system of linear equations 


12°G + (—12n? + 1)C = 12n:@' , 


6n?G —( 12n?+ 1)C =— 6n?@, 
whose solution is 
ve 24n? ., 
(25) C rond , 


(26) ES a & @ 
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where 
(27) De An On ont. 


For the first row (n — 0) we shall find in n. 5 that we must take G = 40), 
where C, is the value arbitrarily chosen for the conductances along the axis. 
Assuming 1/C, = 100000 ohms, we have calculated step by step, using (25) 
and (26), all the resistances of the first 26 rows: they are collected in table I 
along with those for the diagonal network. 

For high values of n we may neglect the last two n —4 3935 
terms in (27), and then formulae (25) and (26) re- 
duce to 


2n ’ 2n +1 
a LÉ mare 
that is 
(28) C=~2nk , 
(29) Ge OE TOE à 


where A is an arbitrary constant. These are the 
formulae adopted by DE PACKH and LIEBMANN. 
On the axis, they give 1/C — oc. The difference 
between these values and those calculated from (25) 
and (26) is far from being negligible in the first two 
or three rows, as is shown in fig. 5. In this we have 
indicated the values of the resistances 1/C, 1/@ 
from our formulae and, in parentheses, those from 
de Packh’s approximation, assuming for the arbi- 
trary constant such a value (K = 3.177 C;) as to 
give coincidence when n tends to infinity. 

The error in the network without diagonals can 
be calculated from (24) (reduced to the third and 
fourth order terms) and (12) with H— H'—0. 
It comes out = 


—1 RDV MOV 
~ 48n? + 2n 6 dr: 24 dr 


(30) € 


Comparing this formula with (22) we see that the diagonals not only sup- 
press entirely the third order error, but eliminate practically even the fourth 
order one, by introdueing in it the very small factor «. 


5. — The first row of resistors (# — 0). 


The calculations of the last two sections are not valid for n — 0, that is 
for the first row of resistors near the axis r — 0. We must start again the 
reasoning, taking into account that, if V is regular on the axis (as we shall 
assume) all its derivatives of odd order with respect to r are zero for r — 0. 
This permits expressing V at any point by means of its values on the axis: 
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in fact, denoting these values by f(z) = V,_,(2), we have the development, 
well known in electron optics (1) 


il 


1 
(31) ne 


So, if we take a junction P, on the axis and the five contiguous junct- 
ions P,, P/, P;, P;, Ps (fig. 6), we have, dropping all terms of higher order 
than the fourth, 


4 


7: if 7 241 h TV 
Vi + Vi=2Ve = hf" st 


h? IOP ops 
i, eee Be 1V 
Ve — Vo= — el “ate ca 
D. 
Re See i meant 


7 1 h? # h = m 13 ra flv 
Vi Vo hf + te: 


(the values of f, f, etc., refer, of course, to point P,). 

In the case of a diagonal network, we may impose the condition that the 
error v— V be zero to the fourth order approximation at all five points: in 
fact Kirchhoft’s equation for point P, gives, with the conductances shown 
10H00, 


: Ce ial Fe CIS : 

2 Y er 24 | | 4 LING = + 
(92) (os ele ag) Bf : 
and it is satisfied for any function f(z) if 
(33) G— 2H =4(,, 

(34) 3G — 26H —=—160,. 


So we can choose at will a value for C,, and determine @ and H by means 
of (33) and (34), which give 


34 
(35) C= = OG esc, 
(36) HE 44 i 


In the case of a network without diagonals, that is with H— 0, we cannot 


(*) See for instance V. E. Cossterr: Bleetron Optics (Oxford, 1950), p. 35. 
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satisfy both equations (33) and (34), and must be content with putting equal 
to zero only the second order part 
of the error: then (33) gives 

(37) Cp == BOW 

and it remains an error of the 
fourth order. 


Formulae (35), (36) and (37) 
have been used for starting the 
computations of table I. 


6. — Boundary meshes. 


We consider now the problem 
of imposing the condition V— V, (V, being a given constant) on a contour 
line ¢ which, generally, will not pass through the junctions of the network. 
Consider, for instance, the case of fig. 9, where three of the eight points 
connected to P, fall out of the boundary c. The network is so calculated that 
the potential at P, would take the right value (in our approximation) if P,, 
P,, P, were at the potentials (5) 


ee yh yee 
(38) Hee Vey Ne 2 
| Vs = Vo EE Ph + (Vee + Ver + 27,)n2/2. 


Now, the current received by P, from P,, P,, P,, in this hypothesis, 
would be 


In order that the same current be carried to P, from a point at potential V,, 
this point must be connected to Py by a conductance S such that 


(39) SE VCH Vo) AV) AV Ve) 


We can realise this condition (without suppressing the connections existing 
in the network) in the following way. If 8, calculated from (39), happens 
to be > 0+ G-+ H (as will usually be the case) we put P,, P,, P, at po- 
tential V, by means of a potentiometer, and add, between them and P,, 


(5) We work out the calculations for the case of cylindrical coordinates; the formulae 
for cartesian coordinates will be obtained simply by replacing the subscripts 2, r with 
x, y and letting the coordinate r tend to infinity. 
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a supplementary resistor of conductance 


(40) S'—= S—(C + G+H), 


as is shown in fig. 7. In the case S< C + G + H, we connect (fig. 8) 
P,, P,, P, (short eireuited) to the 


Bee SA point at potential V,, through a re- 
sistor of conductance 8” given by 


c 
: (41) = 2 1 
i SITES TO CONTE 
In this way, in both cases, the total 
a conductance between P, and the 
0 


point of the potentiometer at po- 
RO EN tential V, will be 8. 


à We now proceed to calculate 8 
3 ne by means of (39). Using (38) and 
Fig. 7 putting V,— V, = U, this formula 
becomes 
Ch h Gh > D 
(42) = U (v. ar Va U (v, ae Vis ci Ar 


Hh 


R 
+5 Le + Ve + Pat Ve 2) J | 


The values of V,, V,, ete., in this formula will be so chosen, that for any 
point of the contour ce, (of coordi- 
nates &, 7 with respect to P,) the 
condition 


(43) V.E + Vin + 
(Vase Vet Verh) 2 0. 


be satisfied as nearly as possible. 
Since the quantities at our disposal 
are five, but are related by equat- 
ion (2), we can satisfy exactly this 
condition only at four points of ec, ; 
which can be chosen at will in the Fig. 8. 
vicinity of Pp. 

It seems quite natural to choose as two of these points Q, and Q, (see fig. 9), 
of coordinates € =h,, 7 — 0 and é=0, 7 = h,: this gives the two equations 


[| Vu V2 =U, 


(44) | V,h, at V,.h2/2 er U : 
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Solving for V,, V, and replacing in (42), taking into account (2), we get 


, PR DV 1 
= 0h} Rh) og + a, je 2 1.) er 
1 VA jemi Ves Vn 
(45) + Gh E a (h hy) an | Hh h, LA | (hy hy) oU +h U 1 
(OU EE À 
+ h, (hk — M), Se = 


Of course, if less, or more, than three junctions connected to Po fall outside 
the boundary, the formula will have less or more terms similar to those written 
In the case of cartesian coordinates the formula reduces to 


R ’ Vay 1 Ving 
(46) S — Ch hi — (h— h,) a + Gh ie + (h—h,) 2 En 
il | V Va 
| h É 2 | ) = uv (2 , Bee 
1 H % + ig (mh) op +h 2 
and since 0 — G, H = (/4 
S 5 Ch + 1 (I h Eh Va 
tr UE = = tv, — My = ee 
(47) A 4 v hy he + 1 se 9 U BT 


The cases of networks without diagonals can be obtained, of course, from (45) 
and (46) by putting H—0. 
We can now proceed to different kinds of approximations. 


a) First approximation. — Neglecting altogether the second order terms 
n (45), this gives 
ei) Ne ler 
ee + x( . 2 
Ve = poem 
hy hy hy, : 


and for the cartesian coordinates, 
with diagonals 


5 ( 3 


and without diagonals 


: h h 
my sehe) 
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Last formula is that indicated by LIEBMANN, and employed also in the re- 
laxation technique. 

This approximation is equivalent to replacing the curve ¢ with a straight 
line between Q, and Q,: it will be sufficient in most cases, if ¢ is not too sharply 
bent in the vicinity of P,. But since the networks with diagonals allow using 
a larger value of h, ir many cases À will be comparable with the radius of 
curvature of c, and a better approximation will be necessary. 


b) Second approximation. — In order to determine V,, and V,, in (45), 
we heed two more conditions, which can be chosen in several ways. One is 
to impose the condition V— V, to point Q, (where ce crosses the bisectrix of 
the axes, fig. 9), of coordinates € = 7 — h;, and to an infinitesimal element 
of the tangent to ce at Q,. Calling A — (dn/d&)o, the slope of this tangent, 
we have, by differentiation of (43), 


1% + AV, ar (We; SF AVer)hs vee À)h; == Or 
and from the condition at Q,: 
We Sr Ve os (1722 Sr Ms SF 2V,r)hs/2 — U/hs 9 


Eliminating V,, V, and V,, by means of (44) and (2), last formulae give the 
system of two equations in V,,, Vz: 


h h, 
(% 2 (1 =| : 


(> hs) Vir = h3(1 ar À) Ver — 


=u 


1 al 1 nef) 
+ (is) rae Fe = 0 | Der |. 


2 


h, — h, 
| > PART ER ART 


Solving for V,, and V,,, and introducing the abbreviation 


A 1 144, A—Ihj—dh, 
ARTE | 

(51) HE MR RER hy 

Me Fs aah (eh ei 
= 


at 


we obtain 


1%; 
52 Ale 
(52) 2 

V, ie el | 1 (=) A | hs 
53 Zu = nee ee N = } si, , 2 
(53) U Ar h m j rh, | his 2 hy + (hs Ms) 5 


These values must be substituted in (45) in order to obtain S and, therefore, 
by (40) and (41), the conductances 8’ or 8” to be put in parallel or in series 
with ©, G, H. In this way, we have taken into account with a nl 
accuracy the shape of the line e between Q, and 0. 

If only two or three junctions are outside the boundary the formula is 
simplified in an obvious manner. 
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In the case of cartesian coordinates with diagonals, the final formula for S 
comes out, from (46), 


5 5 jh il ] 
54 = Y , 
(54) s TOM fe + Ate =) 


(A being given by (51), but without the terms with 1/r) and without diagonals 


5 Ze | 
(55) = A 7} + Au) 


bo 


Other cases of boundary meshes could be treated by the same method. 


7. — Conelusion. 


We believe that a resistor network, constructed as indicated, with the 
diagonal resistors, can give results much more accurate than one of the same 
number of resistors without diagonals, especially in regions where the high 
derivatives of V are appreciable, that is (in the case of electrostatic problems) 
where the field components are rapidly changing. 

We think also that the ordinary network (without diagonals) for cylindrical 
coordinates, will give better results in the vicinity of the.axis if calculated as 
indicated in section 4 rather than with de Packh’s formula. 

In both cases, the method indicated in section 6 enables to introduce the 
boundary condition with greater accuracy, especially in regions where the 
radius of curvature of the boundary is not very large in comparison with the 
mesh side. 

A network with diagonals for cylindrical coordinates is now being con- 
structed in this laboratory. 


The Author is glad to aknowledge the assistance of Miss F. MAGISTRELLT 
in performing the numerical calculations. 


RIASSUNTO 


In questo lavoro si studia la teoria di un tipo perfezionato di rete di resistenze 
per il calcolo numerico di una funzione armonica (in due dimensioni, o in tre con sim- 
metria assiale) con dati valori al contorno. Il vantaggio dell’apparecchio qui studiato 
su quelli finora descritti è che esso tiene conto di termini di un ordine più elevato 
in h (h = lato delle maglie della rete) e permette quindi una precisione molto mag- 
giore con lo stesso valore di h. Si dimostra inoltre che (anche per la rete del tipo usuale) 
il ragionamento basato sulla analogia con la vasca elettrolitica & fallace, e che il metodo 
corretto di calcolo conduce, nel caso della simmetria assiale, a valori delle resistenze 
che, in vicinanza dell’asse, differiscono sensibilmente da quelli adottati da Autori pre- 
cedenti. Si indica infine un metodo per imporre la condizione al contorno sopra una 
linea curva non passante per i nodi della rete. 
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Ultrasonie Grating remaining after stopping the Supersonie Waves 


A. CARRELLI e F. PORRECA 


Istituto di Fisica dell’ Universita - Napoli 


(ricevuto il 15 dicembre 1951) 


During other optical researches on 
the solutions diffusing the light, the 
observations on a colloidal solution of 
starch in water (rice starch), ¢ = 0,1 g/l, 
were showing a peculiar effect which we 
now describe. This solution was put under 
the action of standing supersonic waves 
produced by a quartz, 2 = 0,81 mm; 
great care has been taken to obtain good 
conditions of standing supersonic waves; 
opposite to the quartz was disposed 
a brass blade which could be opportunely 
moved with micrometric movements. 

We measured at first the distribution 
of the intensity of the fringes diffracted 
by the supersonic grating in pure water. 

Table I experimental and 
theoretical values, which are in a sa- 
tisfying agreement: we can think of 
having obtained stationary conditions. 
Putting instead of water the colloidal 
solution of starch in water, we could 
observe that when the circuit generating 
ultrasonics was interrupted, only the 
fringes of odd order disappeared, while 
those of even order remained clearly vi- 
sible for a while (about 6 seconds), gra- 
dually reducing in intensity. 

By means of a phosporoscope and 
using the photographic method, we mea- 
sured the proceeding of the intensity of 


shows 


the first and second order, that remains 
when the ultrasonics stop, at the chang- 
ing of time. 


ABB de 


1 


Order of | Theoretical Experimental | 
the fringes values values 
BT = 3 
0 | 100 100 
1 | 47.5 43.5 + 5.0 
2 | 15.5 17:37 20825 
3 | 6.0 5.8 + 2.0 
4 1.5 1.4 + 1.0 


Reporting the time on the abscisses 
(the zero corresponds to the instant of 
opening of the circuit generator of super- 
sonic waves) and reporting on the ordi- 
nates the intensity, we obtained the 
curves represented in fig. 1. 

First we must observe that this effect 
is strictly connected with the presence of 
standing supersonic waves; aS soon as 
the stationarity is destroyed, or only less- 
ened by the displacing of the brass 
blade, the fringes are not preserved at 
the stopping of the ultrasonics. 


ULTRASONIC GRATING REMAINING AFTER STOPPING THE SUPERSONIC WAVES 91 


It must be remarked that the measur- the considered solution is certainly in- 
ed variation of refractive index An of  ferior to 1076. | 

The distribution of the intensity of 
the fringes in the solution, when the 
ultrasonics are excited, is given by 
table II, which appears different from 


104 


TABLE II. 


Order Experimental 


of the fringes | values | 
0 100 
1 | 52.0 + 5 | 
2 | 25.3 + 3.5 
3 17 32 
4 48 +1 | 


the values above related in tab. I. The 
same effects are present in solution of 
starch in alcool 95%, ¢ = 0,1 g/l; but 
the intensity of the visible fringes when 
the ultrasonies stop, disappears more 
1 swiftly, as may be seen in fig. 1. Those 
effects do not appear in solutions of 
@ starch in benzol. 
All solutions examined are absolutely 


acest invariable throughout the whole expe- 
a) Order 2 Alcool and starch c= 0,1 ; 
db) Order 1 Alcool and starch 
c) Order 2 Water and starch c= 0,1 
d) Order 1 Water and starch c= 0,1 


04” 0,8” 42” € 


/ rience. 
x Further details will be published in 
A a short time. 
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A Probable Example of Multiple Production of Mesons. 


À 


M. DEMEUR, C. 


DitwortH and M. SCHÖNBERG 


Centre de Physique Nucléaire de V Université Libre de Bruxelles 


(ricevuto il 17 Dicembre 1951) 


The star which we here describe was 
found in a group of plates exposed in a 
magnetic field of 22000 gauss at the Pic 
du Midi (2800 m). The emulsion was 
600 micron thick, the Ilford G5 twice 
diluted (1) in order to compensate for the 
low field (?), and was processed by our 
standard methods (?). This star pre- 
sents several interesting features which, 
in our opinion, make it worthy a sepa- 
rate publication. A general outline will 
be given in this short letter and a more 
complete discussion will be published 
later. 

The star (fig. 1) has six minimum 
ionisation branches of small angular di- 


vergence, no visible primary, and à 
single low energy proton track. There 


is a small angle pair 500 micron from the 
centre of the star, nearly parallel to, and 
two micron from, one of the branches. 

The multiple scattering of the tracks 


() ©. WALLER: Bristol Photographie Con- 
gress, 1950. 

(2) C. DILWORTH, 8. GOLDSACK, Y. GOLD- 
SCHMIDT and F. LEvY: Phyl. Mag., 7, 11 (1950). 

() ©. DILWORTH, G. P. S. OCCHIALINI and L. 
VERMAESEN: Bull. Cen. Phys. Nucl. Bruxelles, 
no. 13a; ©. DILWORTH, G. P. S. OCCHIALINI and 
A. Bonerir: Bull. Cen. Phys. Nucl. Bruxelles. 
MO LED. 


was measured by the angular method (4), 
the movement of the stage being con- 
trolled throughout, the distortion of the 
emulsion was measured (5) and found to 
be negligible compared to the scattering 
of the tracks. The sign of the charge 
of five of the branches and of the two 
tracks of the pair could be determined 
with a reasonable probability (2), the 
sixth secondary had too short a track 
length. Grain counting on the six se- 
condaries and the two tracks of the pair 
gave the same value within the errors 
of measurement. The total kinetic energy 
of the secondaries is 9- 10% eV. The 
results are given in Table I. 

The probability of the sign has been 
corrected for distortion, the correction 
in all cases being small. The angle à is 
given with respect to the vertical during 
exposure. 

It seems reasonable to assume that 
this star represents a multiple production, 
rather than plural, given the lack of an 
evaporation star. If we assume the sixth 
unknown charge to be negative, we can 


(4) Y. GoLDSCHMIDT: Bull. Cen. Phys. Nucl. 
Bruxelles. no. 14; Nuovo Cimento. 7, 331 (1950). 

(5) M. Cosyns and G. VANDERAEGHE: Bull. 
Cen. Phys. Nucl. Bruxelles, no. 15. 
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have the conservation of charge required 
for multiple production. We therefore 


40041 


21 de 
Sr 
Fig. 1. — Scale, 5mm=10 u 
P is the low energy proton, a,b, c,d, e, f are 
the fast secondaries, 1 and 2 the two tracks of 
the pair. 


analysed the star on the basis of a 
nucleon-nueleon collision, though we can- 
not be certain that the primary was a 
neutron. 

At first sight our star looks rather like 
a low energy version of the star of 
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SCHEIN (5) but the analysis shows that 
there are important differences. There 
is one evident fact which is diffieult to 
fit into either of the current models of 
meson production (75%). Given thenu m- 
ber of secondaries, their total energy is 
too low for the Fermi model and too 
high for the Heisenberg-Wataghin model. 
Also with such a multiplicity Fermi’s 
model predicts that two thirds of the 
charged particles should be protons and 
antiprotons, whereas no negative proton 
has been detected in this star and all the 
particles could be eventually mesons. 

On the other hand, if we assume that 
all the secondaries are mesons and take 
the energy of the primary to be 6 : 10!!eV 
as required by Fermi’s theory for such 
a multiplicity, we find in the baricentric 
system all the charged particles emitted 
backwards. If we take the primary 
energy to be 1.5- 10! eV as given by 
HEISENBERG for this multiplieity, we find 
that all the charged particles go forward 
in the barycentric Also on 
Heisenberg’s theory the average angle 
would be about 23°, much larger than 
we have. 

In order to have the same number 
of charged particles emitted forwards as 
backwards in the barycentric system, 
we must take for the energy of the pri- 
mary ~7:-10% eV. If the slow proton 
is a recoil, it would then be emitted in 
the barycentrie system with an energy 
of 6: 109 eV, whereas the average energy 
of the mesons would be about 200 MeV. 
There should be an unobserved fast 
neutron with about 5: 10% eV. Accord- 
ing to this picture, only a small part of 
the energy of the nucleons in the bary- 


system. 


(6) J, Lorp, J. FAINBERG and M. SCHEIN: 
Phys. Rev., 80, 970 (1950). 


(7) E. Fermi: Prog. Theo. Phys., 5, 570 
(1950). 

(8) W. HeisenperG: Zeits, j. Phys., 113. 
61 (1939); Nature 164, 569 (1949); VATER 
Phys. 126, 569 (1949). 


() G. WaATAGHIN: Ann Acad. Brasil. Se. 
8, 129 (1943); Phys. Rev, 74, 975 (1948). 
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centrie system would be transferred to 
the meson field. This would seem un- 
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if this star represents a normal or ax 
exceptional case of multiple production, 


Direction Angle of 
Track| #2 plane Dip {scattering 
of plate 0/100 y 
Q p 10-2 
ee ee: 
Bien) p18 15a) Felsen US 
AN ee =e Hagniinaars 
Bo AU ERP 
+206 |—"!o12’| 17 
f | +2018 | +146 | 3,2 
Pair 
à + 2021" | — Me 42% | 713 
2 | +29 |— 1354| 3,0 


likely either on the Heisenberg- Wataghin 
or the Fermi model. 

The mesons are emitted in the bary- 
centric system at rather large angles to 
the direction of the incoming particle. 
With a primary energy of the order of 
7.1010 eV no meson is emitted at an 
angle of less than 45° nor greater than 
135°. This angular distribution is not 
quite what should be expected according 
to either of the two models. It is quite 
possible that the primary was not a 
neutron. If this is so the application 
of models of multiple production due to 
nucleon-nucleon collision may be mis- 
leading, especially if the primary had a 
weak interaction with nucleons, as is the 
case with photons. In order to know 


Grain Energy 

density | Nature |(for Meson) | Sign Proba- 
/100 u 10° eV bility 
30 + 1 <P | 0,80 +1) — LOG 
30 +1 <P | 0,70 +1 | — 7598 
30 +1 ? 1,0 + 2 + 96% 
30 + 1 9 1,7223 + 97% 
32 + 1 ? 1,4 +1 + | 88% 
33 +1 IR, 0,7 + 2 2 — 

for electron 

See e 20.220 | oe eager 
si] e 0,9 +1 — 130% 


it would be desirable to have more 
examples of this type of produetion. 


Acknowledgements 


We wish to thank Professors CosyNns 
and ÖCCHIALINI for their generous as- 
sistance and many useful discussions. 
We are grateful to Prof. GÉHÉNIAU for 
his interest in this work. 

We acknoledge gratefully the assi- 
stance of Mr. WALLER of Ilford Ltd. and 
of the personnel of the Osbervatoire of 
the Pie du Midi, and our colleagues of 
Manchester University, who all contri- 
buted to the task of providing, trans- 
porting and exposing the plates. 


PRES EEC BNZUERTEEPRECENSIONI 


R. M. BozoRTH - Ferromagnetism. 
Un vol. di xvit-+ 968 pp., con 
825 fige., 91 tavole e 4 appen- 
dici. Van Nostrand, New York, 
1951. 


Un trattato sul ferromagnetismo, dove 
i numerosi e complessi aspetti scientifici 
e teenici di questa importante branca 
della Fisica siano trattati abbastanza 
esaurientemente, quindi non a carattere 
monografico, desta sempre interesse: 
questo perchè pochissime sono le opere 
che, trattando esclusivamente del ferro- 
magnetismo, soddisfano a questi re- 
quisiti. 

Merita pertanto d’essere segnalato il 
bel volume che, preparato dal BozorTH, 
che & uno dei fisici più competenti di 
ferromagnetismo e layora con la Bell 
Telephone Co., viene ora pubblicato sotto 
gli auspiei di questo grande complesso 
industriale. 

L’opera risulta divisa in quattro 
parti. La prima, in due capitoli, ha scopo 
introduttivo ed & utile per le definizioni 
fondamentali e per introdurre i non spe- 
cialisti al linguaggio ormai generalmente 
fissato dalla consuetudine in questi studi. 
Scopo puramente informativo e generico 
ha del resto anche la quarta parte, pure 
molto breve, dove sono descritti nel loro 
principio di funzionamento i fondamen- 
tali metodi di misura delle grandezze più 
interessanti per il ferromagnetismo. 

Ben diverso il carattere delle parti II 
e III, di gran lunga più interessanti. La 
seconda parte, in 373 pagine tratta delle 
proprietà magnetiche dei materiali, quin- 
di del Fe, del Ni e del Co, e delle leghe 
di questi metalli tra loro e con altri ele- 
menti, quali il Si, il Mo, il Cr, il Cu, ece.. 


Qualche pagina è pure dedicata alla leghe 
di Heusler. Un intero capitolo è infine 
dedicato ai magneti permanenti e alle 
varie leghe adatte alla loro preparazione, 
come pure alle speciali tecniche che per 
questo scopo sono state elaborate pit 
recentemente (v., ad esempio, metal- 
lurgia delle polveri—magneti preparati 
con polveri di metalli o di ossidi magne- 
tici, pressate, ecc.). 

E questa una parte di notevole inte- 
resse, specialmente per i tecnici e i co- 
struttori, perché in essa si trovano ben 
classificati tutti i tipi di materiali ferro- 
magnetici, con le caratteristiche piu di- 
verse: una grande quantita di utili dia- 
grammi ne indica il comportamento e ne 
fornisce, insieme con le numerose tabelle, 
le principali costanti fisiche, magnetiche 
in particolare. 

La terza parte (di oltre 400 pagine) 
è invece la più interessante per la fisica 
del ferromagnetismo. Dopo aver ricor- 
dato in essa le principali teorie svilup- 
pate per la interpretazione del: ferro- 
magnetismo & per la rappresentazione 
della sua dipendenza dal campo e dalla 
temperatura, un intero capitolo è dedi- 
cato alla teoria dei domini e alla inter- 
pretazione con essa delle particolarita 
della curva di isteresi. Un altro capitolo 
si occupa delle proprietà magnetiche dei 
cristalli. 

Più di cento pagine sono quindi dedi- 
cate alla descrizione dei vari fenomeni 
nei quali più evidenti risultano i feno- 
meni tra magnetizzazione, tensioni mec- 
caniche e deformazioni. In particolare 
sono esaminati i fenomeni di magneto- 
strizione, sia dal punto di vista teorico 
come da quello dei risultati sperimentali : 
ben trattata la questione delle relazioni 
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tra il modulo di elastieita e lo stato di 
magnetizzazione. 

Seguono due capitol in cui sono stu- 
diate le relazioni tra temperatura e ma- 
onetizzazione nei ferromagnetici, con par- 
ticolare riguardo ai fenomeni di carat- 
tere energetico che accompagnano il pas- 
saggio attraverso il punto di Curie. Un 
altro capitolo esamina i fenomeni della 
magnetoresistenza nei cristalli e nei ma- 
teriali policristallini, in relazione anche 
con la temperatura. 

Non molto esauriente appare invece 
il capitolo che si occupa dei cambiamenti 
di magnetizzazione col tempo, in quanto 
non vi sono considerati tutti i tipi di 
«after effects » che sono stati presi in 
esame anche nei recenti lavori di SNOEK. 
Vi è pero considerata la risonanza ferro- 
magnetica. 

Troppo breve e frammentario appare 
infine il Capitolo 18 in cui, dopo aver 
preso in breve considerazione i vari tipi 
di energia magnetica che sono associati 
con la struttura a domini, viene indicata 
la possibilità di prevedere, in base al 
reciproco gioco di queste energie, da 
quali fattori possono essere influenzati la 
permeabilitä iniziale e la forza coercitiva, 

La vasta materia oggetto dell’opera, 
esposta in un inglese corretto e chiaro, 
è svolta in forma molto piana, con assai 
limitato uso di sviluppi matematiei: le 
numerosissime figure, spesso costituenti 
schemi originali, sono disegnate con 
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grande cura e chiarezza. Oltre 2400 lavori 
sono infine eitati per una estensione degli 
argomenti trattati nel libro. 

L’opera del Bozortx, pubblicata 12 
anni dopo il volume Ferromagnetismus 
di BECKER e DÖRING, appare pertanto 
interessante perché aggiornata, e perché 
contiene i dati tecnologici dei diversi ma- 
teriali ferromagnetici, dati che altrimenti 
si dovrebbero rintracciare in libri come 
quello di MmssKin e KussMANN (Die 
ferromagnetischen Legierungen) d’altronde 
vecchio di quasi venti anni. 

Tuttavia, pur nella vastita degli argo- 
menti trattati, il libro di BOZoRTH de- 
lude forse un po’ il fisico, che è portato 
piuttosto a considerare quest’opera come 
una enciclopedia sul ferromagnetismo. 

Devesi perd riconoscere che un mag- 
giore approfondimento della trattazione 
delle singole questioni che interessano il 
ferromagnetismo avrebbe richiesto una 
mole assai maggiore di quella, già rile- 
vante, dell’attuale volume di BozoRTH. 
D’altra parte la trattazione piuttosto sin- 
tetica di una materia vasta e difficile 
come il ferromagnetismo pud presentare 
difficoltà anche maggiori: il modo come 
esse sono state, in generale, superate, 
attesta, se pur ve ne fosse bisogno, della 
competenza che il BOZORTH si è formata 
nello studio del ferromagnetismo. 
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